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The  Naval  Surface  Weapons  Ceriter  Aeroprediction  Code  has  been  extensively 
applied  to  the  prediction  of  static  arid  dynamic  aerodynamics  of  missile  con¬ 
figurations.  Major  extensions  have  recently  been  made  to  the  Code,  extending 
its  capability  to  0  <  M  <  8  and  0~--  ^  a.  <  180d.  and  improving  the  transonic  and 
dynamic  derivative  predictions.  The  theoretical  basis  for  the  Code  is  reviewed. 
The  Code  is  evaluated  through  comparisons  of ^computational  examples  with  experi¬ 
ment  for  body  alone,  body-tail,  and  body-tail-cnnard  configurations.  ,• 
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The  speed  and  accuracy  of  the  Code  is  ideal  for  use  in  preliminary 
design.  Representative  design  charts  (generated  by  the  use  of  the  Code) 
for  both  the  static  and  dynamic  aerodynamic  coefficients  of  a  wide  variety 
of  configuration  components  are  included.  The  charts  will  aid  the  designer 
in  making  preliminary  design  estimates  or  for  data  comparisons. 
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FOREWORD 

This  work  was  undertaken  to  provide  a  design  tool  for  use  in  estimating  the 
aerodynamics  of  tactical  weapons.  Rapid  estimation  of  these  aerodynamics  will 
then  allow  one  to  predict  performance  or  conduct  static/dynamic  stability  analy¬ 
sis  in  the  preliminary  and  intermediate  design  stages  in  an  efficient  and  reason¬ 
ably  acccurate  marner.  Prior  to  the  development  of  this  technology  over  the  past 
10  years,  aerodynamics  were  calculated  by  hand  using  design  charts,  wind  tunnel 
data,  empirical  methods,  or  just  ’’engineering  experience."  It  is  hoped  that  this 
methodology  has  helped  to  lend  an  improved  process  for  generating  aerodynamics. 

Support  for  the  work  war  provided  by  the  following  sponsors:  The  Naval  Sea 
Systems  Command  under  the  Surface-Launched  Weapons  Aerodynamics  and  Structures 
Block  Program,  the  Naval  Air  Systems  Command  under  the  Strike  Warfare  Weaponry/ 
Aerodynamics/Structures  Technology  Block,  the  U.S.  Army  Missile  Command  under 
Project  No.  1L162303A214 ,  and  the  Office  of  Naval  Research  under  Project  No. 
P41411. 

Major  procurement  contracts'  from  the  above  funds  were  let  to  North  Carolina 
State  University,  Nielsen  Engineering  and  Research,  Inc.,  and  Lockheed  Missiles 
and  Space  Co.,  Inc.  Minor  contracts  were  let  to  General  Electric,  Armament  Sys¬ 
tem  Dept.,  Burlington,  Vt.;  C0MPR0  of  Fredericksburg,  Va.;  and  EG&G  of  Rockville, 
Md. 


This  report  was  reviewed  and  approved  by  Mr.  C.  A.  Fisher,  Head,  Weapon 
Dynamics  Division  (K20) . 
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1.  INTRODUCTION 


In  the  design  of  any  airframe,  there  is  a  need  to  estimate  the  aerodynamic 
characteristics  of  the  body  in  flight.  During  the  preliminary  or  conceptual  de¬ 
sign  phase,  design  changes  necessitate  the  revision  of  the  preliminary  estimates 
each  time  a  change  is  made.  To  meet  this  need,  the  Navy  (Army  support  in  1977) 

I 

undertook  in  1971  the  development  of  a  computer  code  which  would  be  able  to  pre¬ 
dict  the  aerodynamic  coefficients  of  a  wide  variety  of  configurations  rapidly, 
inexpensively,  and  with  relatively  good  accuracy.  The  code  was  developed  in 
phases  with  preliminary  versions  of  the  code  being  published  in  1972,  1974,  and 
1977  (see  References  ,  2,  and  3,  respectively).  These  versions  were  for  the 
body  alone  static  aerodynamics,  wing-body-tail  static  aerodynamics,  and  dynamic 
derivatives,  respectively.  Applicable  Mach  number  and  angle-of-attack  ranges 
were  0  <  M^  <  2.5  and  0  <  a  <  stall,  respectively.  The  final  version  of  this 
code,  simply  referred  to  as  the  Aeroprediction  Code,  was  completed  in  late  1981. 
This  version  of  the  code  extended  the  Mach  number  range  to  6  and  aagl.e-of-attack 
to  180°.  However,  wing-body  aerodynamics  can  only  be  computed  to  a  =  45°.  Above 
a  -  45°,  body  alone  aerodynamics  are  generated.  Theoretical  documentation  of 
this  last  version  of  the  code  is  given  in  Reference  4.  The  code  has  been  devel¬ 
oped  to  handle  a  general  class  of  axisymmetric,  wing-body-tail  configurations 
(this  does  not  include  inlets  or  base  flow).  The  general  approach  in  the  code 
development  has  beer,  to  combine  existing  and  newly  developed  computational  meth¬ 
ods  into  a  single  computer  program.  The  basic  method,  described  in  more  detail 
in  Section  2,  is  that  of  component  superposition  where  the  body  alone,  lifting 
surface  alone,  and  interference  contributions  are  added  to  obtain  the  total  con¬ 
figuration  aerodynamics . 

The  purpose  of  this  repute  is  twofold:  1)  to  provide  a  brief  summary  of  the 
methods  used  in  the  Code,  end  2)  to  provide  representative  design'  charts  which 
have  been  generated  by  the  use  of  tne  Aeroprediction  Code  and  which  can  be  used 
for  preliminary  design  estimates  or  data  comparison.  Section  2  provides  the  sum¬ 
mary  cF  the  prediction  methods.  The  methods  are  outlined  in  sufficient  detail  so 


that  one  familiar  with  aerodynamic  theory  will  probably  not  need  to  seek  further 
detail  in  the  earlier  works.  However,  many  of  the  intermediate  steps  and  data 
tables  used  in  the  Code  are  not  presented.  Section  3  provides  a  comparison  of 
the  theory  in  the  Code  and  experimental  data  for  a  number  of  representative 
configurations  which  will  provide  the  users  of  both  the  Code  and  the  design 
charts  a  feel  for  the  overall  accuracy  of  the  predictions.  There  are  also  sec¬ 
tions  discussing  the  design  of  projectiles  and  missiles  to  achieve  optimum  per¬ 
formance.  Section  4  presents  a  limited  number  of  design  charts  which  can  be  used 
by  the  designer  to  obtain  an  estimate  of  the  aerodynamic  coefficients  for  a 
variety  of  configuration  components.  The  chief  purpose  of  the  charts  is  to  show 
the  trends  in  the  coefficients  which  result  from  changes  in  a  configuration 
component.  The  charts  in  many  cases  represent  reproductions  or  extensions  of 
charts  provided  in  earlier  works ,  which  are  shown  both  for  comparison  and  to 
provide  new  data.  The  charts  presented  by  no  means  represent  the  full  range  of 
application  of  the  Aeroprediction  Code.  Time  and  funds  would  not  allow  a  com¬ 
plete  examination  of  the  effects  of  varying  geometric  parameters.  The  user  may 
wish  to  conduct  additional  parametric  analyses  within  the  applicable  limits  of 
the  Code. 


2.  METHODS  OF  ANALYSIS 


The  basic  approach  is  that  of  component  linear  superposition.  Body  alone, 
lifting  surface  alone,  and  interference  aerodynamic  contributions  are  combined 
into  total  configuration  aerodynamics.  No  inlet  or  plume  effects  are  considered. 

The  general  approach  has  been  to  combine  existing  and  newly  developed  ap¬ 
proximate  computational  methods  into  a  single  program  to  compute  aerodynamics. 
Computational  times  required  for  the  estimate  of  static  and  dynamic  aerodynamic 
coefficients  for  a  body-tail-canard  configuration,  which  is  the  most  complex 
configuration  considered,  are  in  Central  Processing  Unit  (CPU)  seconds  on  the 
CDC  6700  computer  as  opposed  to  minutes  or  hours  required  for  more  detailed 
physical  and  numerical  models.  The  accuracy  obtained,  however,  is  compatible 
with  that  required  for  preliminary  or  intermediate  design  estimates. 

2.1  STATIC  AERODYNAMICS 

2.1.1  BODY  ALONE  AERODYNAMICS 

Complete  generality  of  body  configuration  is  not  possible.  The  basic  ge¬ 
ometry  for  the  most  complex  configuration  considered  is  shown  in  Figure  1.  The 
body  is  axisymmetric. 
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Figure  1.  General  Body  Geometry 
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The  nose  may  he  pointed,  spherically  blunted,  or  truncated.  For  the 
spherically  blunted  case,  the  origin  of  coordinates  is  where  the  spherical  cap 
joins  the  rest  of  the  nose.  Note  that  x  is  a  coordinate  from  the  nose  tip.  The 
nose  may  also  consist  of  two  piecewise  continuous  sections  in  addition  to  a  pos¬ 
sible  spherical  cap.  At  high  supersonic  riach  numbevs,  the  spherical  cap  may  be 
generalised  to  a  conic  section  curve.  The  afterbody  consists  of  a  constant 
radius  section.  The  boattail  or  flare  may  be  a  conical  frustrum  or  a  curved 
section.  A  rotating  band  for  a  projectile  can  also  be  included. 

The  reference  area  is  associated  with  a  reference  body  diameter  for  all 
aerodynamic  coefficients.  The  diameter  must  be  that  at  the  end  of  the  nose. 

2. 1.1.1  AXIAL  FORCE 

Axial  force  components  are  broken  down  into  linear  superposition  of  wave  or 
pressure  force,  friction  force,  rotating  band  force,  and  base  force. 

2. 1.1. 1.1  WAVE  OR  PRESSURE  AXIAL  FORCE 

The  wave  or  pressure  axial  force  is  defined  as  the  forebody  axial  force  due 
to  normal  stress  or  pressure.  In  the  subsonic-transonic  flow  region,  the  fore¬ 
body  surface  pressures  are  affected  by  boundary  layer  separation  at  the  nose- 
afterbody  junctions,  at  the  end  of  the  body,  and  on  a  steep  boattail. 

2. 1.1. 1.1.1  SUBSONIC  AND  TRANSONIC  FLOW  (Mw  <  1.2).  Currently  the  wave  axial 
turce  is  assumed  to  be  contributed  by  the  nose  and  the  boattail  in  a  linear 
superposition.  For  bodies  with  long  afterbodies,  the  nose  and  boattail  drag  are 
essentially  uncoupled.  In  other  words,  the  nose  axial  force  does  not  depend  on 
boattail  parameters  and  afterbody  length  (subsonic  Mach  numbers),  and  the  boat- 
tail  axial  force  does  not  depend  on  nose  parameters  and  afterbody  length.  In 
reality,  the  boattail  axial  force  and  base  force  are  coupled  for  subsonic  Mach 
numbers,  even  for  long  afterbodies. 

Complete  generality  of  solution  would  require  a  solution  of  small  dis¬ 
turbance  potential,  fu^l  potential,  or  Euler  model  equations.  Computational 
times  are  prohibitive  even  for  the  simplest  small  disturbance  model. 


For  the  nose  axial  force,  the  approach  is  to  utilize  computational  axial 
force  coefficient  data  and  experimental  data  for  a  nose  followed  by  a  long  after¬ 
body  in  a  table  lookup  manner.  Reference  5  indicates  that  the  nose  axial  force 
coefficient  is  shape  dependent.  The  family  of  tangent  ogives  and  conical  noses 
approximate  the  boundaries  of  the  envelope  of  shapes  of  interest. 

The  nose  axial  pressure  or  wave  drag  force  coei  ficient  is  given  by  a  linear 
interpolation  between  a  blunted  ogive  value  and  a  blunted  cone  value. 
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AN  (Rg  +  Rq) 

Details  and  evaluation  of  the  prediction  are  given  in  Reference  4. 

CANO  ''■S  t*ie  tan8ent  ogive  value  and  C^q  is  the  shar  .e  value.  Rq  is  the 
slope  at  x  =  0  and  Rg  is  the  slope  at  the  nose-afterbody  uiction.  When  Rg  is 
zero,  the  ogive  value  is  obtained,  and,  when  Rg  =  Rq,  the  cone  value  is  obtained. 
£g  is  the  cone  angle  associated  with  Rg. 

The  value  for  (R^,  Ljj,  M^)  is  obtained  from  a  table  lookup  interpola¬ 

tion  routine  for  blunted  tangent  ogive  noses  with  long  after-bodies.  Tabular 
values  are  obtained  from  solutions  of  the  full  potential  equation  and  a  model  of 
the  Euler  equations.  The  original  computational  method  reported  in  Reference  6 
obtains  a  solution  of  the  Euler  equations  model  by  an  implicit  time  asymptotic 
method.  Computational  accuracy  from  Reference  6  for  shorter  nose  lengths  was 
deemed  to  be  insufficient.  Hence,  the  full  potential  model  of  Reference  7  was 
used  to  update  the  table.  The  potential  model  data  transmitted  by  G.  Kuhn  cf 
NEAR,  Inc.  are  unpublished.  Discrete  values  of  are  .8,  .95,  1.05,  and  1.2. 
R^j  values  are  spaced  in  .05  increments  from  0  to  .5.  values  range  from  .75  to 
5.  For  values  below  =  .8,  C^q  is  decayed  quadratically  from  the  value  at 

tl  =  .8  to  0  at  II  =  .  5 . 

00  00 

The  tabular  values  were  obtained  by  integrating  pressure  data  for 

pointed  cone-cylinders  from  Reference  8  and  blending  these  data  with  values  from 
the  classical  Taylor-Maccoll  solutions.  va*ues  are  decayed  to  a  finite 

value  asymptote  below  =  .5.  Reference  9  indicates  that  the  effect  of  blunting 
is  not  significant  for  bluntness  as  great  as  R^  =  .4  for  pointed  conical  noses 


below  shock  attachment  Mach  numbers.  £g  tabular  values  are  in  5°  increments 
ranging  from  0°  to  20°.  M^  tabular  values  are  .5,  .7,  .8,  .95,  1.05,  and  1.2. 

The  boattail  axial  force  model  is  reported  in  Reference  1.  For  supersonic 
Mach  numbers  close  to  1.0,  Reference  10  gives  an  approximate  solution  of  the 
small  disturbance  potential  equation  on  a  boattail  which  follows  a  long  after¬ 
body.  For  Mach  numbers  below  1.05,  the  boattail  drag  is  decayed  to  0.0  in  a 
linear  manner  at  M^  =  .95  from  the  value  at  M„  =  1.05. 

oo  oo 

2. 1.1. 1.1. 2  LOW  SUPERSONIC  MACH  NUMBERS  (1. 2  <  M^  <  M£).  M£  is  a  Mach  number 
separating  low  supersonic  inviscid  computations  from  high  supersonic  computations 
for  the  body  alone  components  for  static  coefficient  prediction.  M^  is  a  program 
input  between  2  and  3.5  which  is  dictated  by  the  limits  of  the  small  disturbance 
models  used  for  the  low  supersonic  Mach  number  flow  region. 

Detailed  theory  and  evaluation  for  this  section  are  given  in  Reference  1. 
For  pointed  bodies,  the  basic  model  is  the  Van  Dyke  second-order  perturbation 
solution  of  the  potential  equation  (Reference  11).  The  full  nonlinear  poten¬ 
tial  equation  can  be  separated  into  an  axial  flow  nonlinear  problem  and  a  linear 
nonconstant  coefficient  (constants  depend  on  an  axial  solution)  crossflow  problem 
when  the  angle  of  attack  is  small  (15°  or  less). 

First-order  linear  solutions  are  obtained  by  neglecting  all  nonlinear  terms 
for  the  axial  flow  and  crossflow  problems.  The  first-order  axial  flow  problem  is 
solved  by  superposing  axial  source  solutions  in  integral  form.  The  source 
strength  distribution  is  estimated  by  a  piecewise  second-order  Taylor  series 
representation.  The  classical  KArraAn-Moore  source  distribution  is  first-order 
piecewise  linear.  Van  Dyke  estimates  the  highest  order  nonlinear  terms  using  the 
first-order  solution.  Thus  a  solution  of  the  nonlinear  axial  flow  problem  is 
obtained  by  a  second-order  iteration  or  perturbation.  The  second-order  partial 
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differential  equation  has  the  same  linear  terms  and  is  nonhomogeneous .  An  ap- 
proximate  closed  form  particular  solution  was  found  by  Van  Dyke  to  the  same  order 
of  accuracy  as  the  nonlinear  terms  retained.  The  second-order  homogeneous  solu¬ 
tion  is  solved  using  a  piecewise  seconds  rder  axial  source  distribution. 

To  be  consistent,  the  crossflow  solution  contributes  terms  of  order  sin2  a 
to  the  axial  force  coefficient.  However,  tlie  Code  uses  the  full  nonlinear  pres¬ 
sure  coefficient  relation.  The  crossflow  solution  contributes  a  term  of  order 
sin  a  to  the  normal  force.  The  crossflow  solution  will  be  discussed  in  the  sec¬ 
tion  on  normal  force  cr  lift  at  zero  angle  of  attack.  Only  a  first-order  linear 
crossflow  solution  is  considered  in  the  "hybrid"  model. 

Thus  far,  the  discussion  has  dealt  with  a  pointed  body  solution.  The  origin 
of  the  supersonic  axial  source  and  doublet  distributions  coincide  with  x  a  0. 
The  perturbation  solution  is  limited  by  the  requirement  that  the  Mach  angle  be 
greater  than  the  maximum  body  slope.  In  the  Code,  the  requirement  is  that  the 
maximum  body  slope  be  123S  than  .95  of  the  slope  of  the  Mach  angle  for  pointed 
bodies . 

For  a  spherically  blunted  body  the  slope  exceeds  the  Mach  angle  slope  on 
part  of  the  sphere.  A  modified  Newtonian  pressure  distribution  is  matched  with 
the  perturbation  theory  at  some  point  on  the  sphere.  The  modified  Newtonian 
pressure  distribution  is  given  as 

C  =  C  sin2  6  (2) 

P  P° 

sin  6  =  sin  e  cos  a  -  cos  e  cos  0  sin  a  (3) 

where 

6  is  the  azimuthal  angle  from  the  leeside, 

e  is  the  angle  associated  with  the  body  slope,  and 

CpQ  is  the  stagnation  pressure  coefficient  behind  a  normal  shock. 

The  method  for  matching  pressure  distribution  and  for  selecting  the  origin 
of  the  source  distribution  is  complex.  la  order  to  improve  the  recompres3ion 
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predicted  at  higher  Mach  numbers  in  the  neighborhood  of  the  sphere-body  junc¬ 
tion,  the  nose  shape  was  modified  as  shown  in  Figure  2. 
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Figure  2.  Spherically  Blunted  Nose  Geometry 


The  tangent  to  the  body  at  the  origin  of  coordinates  is  extended,  and  a 
perpendicular  is  dropped  from  the  origin  to  point  t.  For  low  supersonic  flow 
the  spherical  nose  cap  is  replaced  by  a  sphere  and  a  conical  frustrum  recom¬ 
pression  ramp.  The  effective  sphere  radius  is  then  R^  *  RQ  cos  eQ  as  opposed 
to  the  one  shown  in  Figure  1  where  =  R^/cos  e^. 

The  perturbation  solution  starts  at  point  s.  The  point  s  initially  is 
chosen  where  •  the  surface  tangent  angle  is  27.5°.  An  additional  requirement  is 
that 

F  tan  p  <  Rq  (4) 

where 

F  is  an  input,  to  the  code  of  about  .95,  and 
p  is  the  Mach  wave  angle. 

If  Equation  (4)  is  not  satisfied,  then  the  effective  free  stream  Mach  numbv . ..  i„. 
cut  back  to  the  point  where  F  tan  p  =  R^.  If  tan  (27.5°)  <  Rq,  then  s  coincides 
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with  t.  The  origin  for  the  source  or  doublet  distribution,  st,  is  obtained  by 

ew 

extending  a  tangent  from  s  to  the  axis  at  st.  Thus,  the  perturbation  solution 
starts  as  a  cone  solution  and  marches  from  s  down  the  body.  The  modified  New¬ 
tonian  solution  starts  at  the  stagnation  point.  On  seven  meridian  planes  of  con¬ 
stant  8  (0°,  30°,  60°,  ...,  180°),  an  intersection  between  the  two  distributions 
is  sought.  If  no  intersection  is  obtained,  then  a  pressure  discontinuity  is 
allowed  at  s.  However,  even  if  there  is  an  intersection,  there  is  a  pressure 
slope  discontinuity.  The  Newtonian  pressure  distribution  is  used  upstream  of  the 
point  of  intersection.  The  perturbation  solution  is  used  downstream  of  the 
intersection. 

For  the  truncated  body  case  (Figure  1  with  a  flat  nose),  a  pseudo  sphere 
is  used  to  provide  a  solution  at  the  spherical  cap-body  junction.  The  "hybrid" 
solution  begins  at  point  s  and  continues  down  the  body.  Loads  are  not  in¬ 
tegrated  on  the  pseudo  sphere.  No  pressure  intersection  is  required.  On  the 
flat  face,  the  mean  pressure  is  assumed  to  be  .9  of  the  stagnation  value. 

Once  the  Cp  distribution  is  obtained  from  the  "hybrid"  model,  the  axial 
force  coefficient  is  given  by 

L  n 

C*  =  8/n  [  f  C  (x,0)  RR'dQdx  .  (5) 

A  0  \  p 

Dimensions  here  are  in  calibers.  Cp  is  given  by  the  full  nonlinear  isentropic 
Bernoulli  relation. 

A  lower  limit  criterion  for  is  given  approximately  by 

.95  tan  p  >  R' 

r  max 

For  a  spherically  blunted  body,  is  for  the  part  of  the  body  downstream  of 

the  spherical  cap.  Numerical  stability  and  accuracy  is  degraded  for  blunt  bodies 
as  the  Mach  number  increases. 

2. 1.1. 1.1. 3  HIGH  SUPERSONIC  FLOW  (M^  >  M£).  Detailed  theory  and  evaluation  for 
this  section  is  given  in  Reference  A.  The  method  used  is  based  on  a  modification 
of  that  given  in  Reference  12. 
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The  body  in  Figure  1  consists  of  piecewise  continuous  elements.  Ths  first 
element  is  a  conic  section  of  the  form 

R2  =  2  Cx  +  Bx2  .  (6) 


The  remainder  of  the  body  consists  of  two  nose  sections,  an  afterbody,  and 
a  boattail  or  flare.  Some  r.f  these  sections  may  be  omitted.  Each  piecewise 
continuous  section  is  divided  into  n  (currently  40)  increments.  This  defines 
points  where  the  pressure  is  to  be  calculated.  Tangents  are  drawn  to  the  body 
at  these  points  and  extended  until  intersections  with  neighboring  tangents  are 
found.  Thus,  the  body  is  replaced  by  a  circumscribed  body  of  linear  segments. 
On  a  pointed  or  truncated  nose  the  first  tangency  point  is  at  the  origin.  On  a 
conic-section  blunt  nose,  the  first  tangency  point  is  determined  by  a  match  with 
e  modified  Newtonian  pressure  distribution.  If  the  first  section  of  the  nose 
is  a  pointed  cone,  obviously  no  subdivision  is  necessary. 


The  pressure  coefficient  distribution  is  assumed  to  be  represented  as 

C  =  (C  V  „  -  A  sin  2o(  cos  0  +  (F  cos2  0  +  A  sin2  0)  sin  2  a/2  (7) 

p  p  a=0 

where  A,  f,  and  A  can  be  represented  as 


A  =  -  I  (  — E 

2  V  3a 


92C 

r  =l 


a=0,  0=0 

,  and 


a=0,  0=0 


/92C  , 

4  =  l  ss^L 


a=O,0=n/2 


(8) 

(9) 

(10) 


From  second-order  shock-expansion  theory,  which  neglects  incoming  characteristics 
at  the  body  surface  and  makes  other  approximations  for  the  outgoing  characteris¬ 
tic  equations,  one  may  show  that  when  (C  )  -,  the  A,  T,  and  A  variations  on  a 

p  a-u 

straight  line  segment  are  given  by  the  same  functional  form 
Y  =  Ya  +  (Y2  -  Ya)  e‘n  ,  and 
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Here  s  is  the  distance  along  the  linear  segment,  subscript  2  refers  to  the  value 
at  the  beginning  of  the  segment,  and  subscript  a  refers  to  an  asymptotic  value. 

Slopes  are  discontinuous  across  a  segment  junction.  Jump  conditions  across  a 
slope  discontinuity  are  obtained  from  Prandtl-Meyer  isentropic  relations  and 
higher  derivatives  thereof.  Pressure  coefficient  values  just  after  a  corner  can 
thus  be  related  to  values  just  before  a  corner  (indicated  by  the  subscript  1). 
Asymptotic  values  are  given  by  an  approximate  cone  solution  for  positive  values 
of  local  e  and  free  stream  Mach  number.  The  pressure  gradient  jump  at  a  corner 
is  given  by  a  relation  from  Reference  12. 

For  negative  values  of  £,  *s  assume<l  to  be  zero  and  the  remain¬ 
ing  asymptotic  values  of  the  pressure  functions  are  given  by  a  blend  of  slender 
body  and  hypersonic  approximations  at  free  stream  Mach  numbers. 

For  a  blunt  body,  the  initial  values  of  the  p  ’sure  coefficients  are  ob¬ 
tained  by  evaluation  of  Equation  (2)  and  its  derivatives  at  the  match  point. 
The  match  point  is  where  the  local  Mach  number  is  1.1  as  determined  by  a  com¬ 
bination  of  Equation  (2),  evaluated  at  a  =  0,  and  the  isentropic  pressure  rela¬ 
tion.  For  a  pointed  body  or  truncated  body  the  initial  values  of  the  pressure 
functions  are  obtained  from  the  approximate  cone  solution  evaluated  at  £q  at 
the  free  stream  Mach  number.  For  the  truncated  flat  face,  the  average  pressure 
is  again  assumed  to  be  .9  of  the  stagnation  value.  Values  at  the  tangency  point 
and  at  the  end  of  the  segment  are  obtained  from  Equations  (11)  and  (12)  for  the 
pressure  functions. 

Once  the  pressure  distribution  is  determined,  is  obtained  from 
L 


C.  =  8 
A 


{  [(V-0  + 


r+A  .  2 

—  sin 


ajdx 


(13) 
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2. 1.1. 1.2  FRICTION  AXIAL  FORCE 

The~ boundary  layer  on  the  surface  of  an  open-end  cylinder  is  analogous  to 
that  on  a  flat  plate  of  the  same  length.  The  model  used  here  is  that  of  Van 
Driest  (Reference  13)  and  is  given  in  detail  in  Reference  1.  A  portion  of  the 
forward  part  of  the  plate  is  assumed  to  be  laminar.  A  critical  Reynolds  number 
based  on  body  length  is  chosen  as  106  (Reference  1). 

The  friction  axial  force  coefficient  is  given  as 


“AF  "  n 


where 


is  the  body  surface  area  in  calibers,  and 

Rjjk  and  R^c  are  Reynolds  number  and  critical  Reynolds  number,  respectively, 
based  on  body  length. 


For  Rnl  <  Rnc, 
CAF  =  rt  Sb 


1.328 

tFnl 


from  the  well-known  Blasius  flat  plate  solution.  Cj,  and  Cp,^  are  friction  co¬ 
efficients  obtained  at  the  Reynolds  number  and  critical  Reynolds  number,  respec¬ 
tively.  In  general,  Cp,  is  given  by 

di/  ifc^T  -  log  CF  =  d2  (16) 


where 


di  =  .55723  [sin"1  (Cj)  +  sin"1  (C2)]/fT/ T^~, 
d2  =  log  RNL  -  1.26  log  (TJTJ  , 

„  _  2A2  -  B 

'■a - r  t 

(B2  +  4A2)* 


(B2  +  4A2)* 
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w/T 


00 


B  = 


1  +  M» 
_ 2  » 


-  1 


and 


VT*  =  1  +  .9  ^  Hi  . 


00 


Equation  (16)  is  solved  twice  by  a  Newton-Raphson  iteration  algorithm  for 
the  input  and  R^  =  R*^. 

2. 1.1. 1.3  ROTATING  BAND  AXIAL  FORCE  COEFFICIENT 

The  band  height,  H,  is  assumed  to  be  much  less  than  the  band  width.  The 
band  is  also  assumed  to  be  far  back  on  a  constant  cross  section  afterbody  in 
fully  developed  turbulent  flow.  j 


Then 


CARB  =  ACA  (fU  H/'01 


(17) 


where  H  is  the  band  height  in  calibers  and  AC^  (M^)  is  a  table  lookup  axial  force 
coefficient  for  a  band  .01  caliber  in  height. 

2. 1.1. 1.4  BASE  AXIAL  FORCE 

It  is  assumed  that  the  boattail  or  flare  is  located  after  a  relatively  long 
afterbody  so  that  the  approaching  external  flow  is  at  free  stream  conditions. 

The  mean  base  pressure,  then,  is  a  function  of  boattail  length,  shape,  and 
final  diameter.  The  model  adapted  in  Reference  1  is 


CAB  =  _CPAB  (RB/Rref)3 


(18) 


where  Cp  (M  ,  1)  is  the  axial  pressure  coefficient  for  a  cylinder  without  a 

.  «-«•  -1  AB 

boattail. 


For  a  flare,  the  model  used  is 

CAB  =  "CPAB  ^VRref^ 


(19) 
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Here.  P  ,  is  the  constant  radius  of  the  afterbody  and  R„  is  the  base  radius. 

C^g  ;*'s  avera8e  base  force  coefficient  from  experimental  data  for 

a  constant  flare  angle  (0  <  8g  <  15°).  Analysis  of  the  available  data  leads  to 
the  conclusion  that  Equation  (18)  is  approximately  valid  at  higher  supersonic 
Mach  numoers  and  invalid  at  other  Mach  numbers  except  for  very  small  boattail 
angles . 

At  angle-of-attack,  an  incremental  value  is  added.  The  expression  is  based 
on  data  from  Reference  14. 

ACABA  =  (*6493a  ‘  -002833MJ  (RB/Rref)0  (2°) 

Here  a  is  the  angle-of-attack  in  radians  Note  that  the  equation  has  been  re¬ 
vised  from  that  given  in  Reference  X  to  fit  the  data  more  accurately.  For 
flares,  (Rg/Rrefj2  *-s  use^  in  Equation  (20).  Equation  (20)  is  valid  for  the  body 
alone  case.  In  the  section  on  lifting  surface  and  interference  aerodynamics, 
Equation  (20)  is  replaced  by  a  different  expression.  The  total  base  axial  force 
coefficient  is  given  by  CAB  +  ACABA ' 

2. 1.1. 2  NORMAL  FORCE  AND  PITCHING  MOMENT 

The  body  normal  force  is  broken  down  into  inviscid  end  viscous  crossflow 
contributions.  Moment  centers  are  defined  with  respect  to  the  nose.  Reference 
length  is  the  body  reference  diameter.  The  final  moment  is  about  the  moment 
center. 

2. 1.1. 2.1  INVISCID  NORMAL  FORCE 

2. 1.1. 2.1  \  SUBSONIC  AND  TRANS0NIC  FLOW  (hM  <  1.2).  In  Reference  1,  CN  and  Cm 
are  assumed  to  be  contributed  by  a  linear  f "perposition  of  component  par?s  of  the 
body.  Each  component  ’s  assumed  to  be  independent  of  the  other.  The  nose  normal 
force  is  predicted  om  a  fit  of  the  cone  alone  data  of  Reference  ?.  The  center 
of  pressure  prediction  is  given  by  slender  body  theory.  The  afterbody  is 
given  as  a  table  lookup  function  of  afterbody  length  and  free  stream  Mach  number. 
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Both  Cv  and  x  _  are  given  by  the  transonic  small  disturbance  theur,  of  Refer- 
ence  15.  For  low  Mach  numbers,  C„  is  predicted  from  a  limited  data  correlation. 

Of 

For  low  Mach  numbers,  the  x  prediction  is  assumed  to  be  the  same  as  for  tran¬ 
sonic  flow. 


The  boattail  Cjj  is  based  on  limited  data  correlation  where  /[I 
-  (Rg/Rj-gf)2}  is  fittecf  as  a  function  of  "||  |M£-l|  .  x^  is  also  given  by  a  slender 
body  theory  estimate.  The  correlation  functional  form  is  based  on  slender  body 
theory. 

This  combination  of  methods  is  necessitated  by  the  lack  of  a  theoretical 
and  experimental  data  base. 

The  approach  of  Reference  16,  which  was  evaluated  in  some  detail  in  Refer¬ 
ence  4,  was  intended  to  improve  the  predictive  capability  in  transonic  flow 
Complete  generality  obviously  is  not  possible.  The  nose  was  assumed  to  be  a 
blunted  tangent  ogive.  The  boattail  was  assumed  to  be  conical. 

The  individual  body  part  contributions  are  superimposed  linearly.  Howeve”, 
the  individual  contributions  are  assumed  to  be  functions  of  all  body  parameters 
such  as  nose  spherical  radius,  nose  length,  afterbody  length,  and  boattail  angle. 
The  data  base  for  Cjj  and  is  generated  by  solving  the  Euler  equations  at  1° 
angle  of  attack.  Un¥ortunate¥y,  the  size  of  the  data  base  was  not  large  enough 
to  properly  represent  this  limited  model.  Ptrts  of  the  data  bane  wen  also 
suspect  as  to  accuracy,  in  particular  the  boattail  contribution  at  all  Mach  num¬ 
bers  and  the  C„  prediction  at  supersonic  Mach  numbers. 

For  each  individual  part,  cr  C  are  fittod  to  a  quadratic  Taylor  -,.rries 
functional  form  for  four  Mach  numbers  (.75,  .9,  .95,  1.?).  15-fit  coefficients 
are  involved  in  the  Taylor  series  in  four  geometric  p-aameter^.  For  each  Mach 
number,  the  Taylor  series  for  the  nose,  the  afterbody,  and  .5,  1-,  and  1.5  cali¬ 
ber  boattail  lengths  are  used.  The  zero  length  t^zero  contribution),  .5,  1.,  and 
1.5  caliber  estimates  are  used  to  provide  data  for  a  four-point  Lagrange  fit  i or 
the  boattail  contribution.  At  Mw  ~  .6,  the  individual  tody  part  contributions 
are  givtn  by  the  model  from  Reference  1.  A  five-pome  Lagrange  fit  iu  Mw  yields 
the  final  component  part  contributions  in  the  range  .6  <  .1  <  1-2.  For  Mw  <  .6, 
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the  prediction  is  given  by  the  older  model.  The  data  at  =  1.2  have  been  re¬ 
placed  by  a  fit  to  a  large  set  of  data  generated  by  the  supersonic  small  dis¬ 
turbance  Van  Dyke  potential  model  of  Reference  17. 

Currently  either  of  the  two  models  is  an  input  option.  For  afterbody 
length?  of  more  than  5.0  calibers,  the  older  model  is  used  since  the  data  base 
lot  the  more  recent  model  is  limited  to  an  afterbody  of  5.0  calibers. 


2. 1.1. 2. 1.2  LOW  SUPERSONIC  FLOW  (1.2  <  M  <  Mn).  A  line  doublet  distribution 
on  the  axis  p.-rmits  a  solution  of  the  first-order  linear  crossflow  equation.  A 
boundary  layer  thickness  based  on  flat  plate  theory  is  added  to  the  original  body 
before  the  crossflow  equation  it  solved.  As  indicated  in  the  section  on  axial 
force  prediction,  is  based  on  a  '‘hybrid"  combination  of  first-order  crossflow 
velocity  and  second-order  axial  flow  velocity.  The  pitching  moment  about  the 
nose  and  the  normal  force  are  then  given  by 

ft  L  n 

CM  =  -  -  f  f  C  (x,  0)  R  eos&  d0  dx  and,  (21) 

N  n  *o  0  P 

0 

C  =  l  if  x  C  (x,  0)  R  cose  de  dx.  (22) 

m  71  n  \  P 


2. 1.1. 2. 1.3  h  lGh  SUPERSONIC  FLOW  (M^  >  M^).  The  basic  model  has  been  given  in 
the  corresponding  section  on  axial  force  prediction  where  and  Cm  may  be  given 
in  terras  of  the  pressuv  coefficient  as 

L 

Cw  =  4  sin  2a  f  R  A  dx  and  (23) 

N  '0 


C 

m 


-  4  sin  2a 


f  x  R  A  dx  . 
0 


(24) 


2. 1.1. 2. 2  VISCOUS  CROSSFLOW  MODEL 

The  model  of  Allen  and  Perkins  (Reference  l?)  is  used  and  is  given  in  con’ 
den3ed  form  in  Reference  1.  The  theory  is  analogous  to  slender  body  theory  in 
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that  the  normal  force  loading  is  assumed  to  be  given  for  a  local  two-dimensional 
cylinder  exposed  to  the  crossflow  velocity,  Va  sin  a. 

An  incremental  normal  force  coefficient  is  then  given  as 


acn  =  ^r1  2  f  RCdc(M«s'aQ)  d*lsin2a 

•‘0  J 

=  M).  c ,  (Mm  sin  a)  S  sin2  a  . 
dc  »  P 


(25) 


where  C^c  is  an  infinite  cylinder  crossflow  drag  coefficient  given  by  a  table 
lookup  interpolation.  Fully  developed  turbulence  independent  of  crossflow 
Reynolds  number  is  assumed,  r]  is  a  finite  length  correction  factor  which  is 
given  by  a  table  lookup  interpolation,  and  is  the  planform  area  in  calibers. 

Similarly  the  incremental  moment  about  the  nose  is  given  as 

AC  =  -  ^  n  C,  x  sin2  a  (26) 

m  n  1  dc  p 


where  x^  is  the  center  of  the  planform  area. 

2.1.2  LIFTING  SURFACES  AND  INTERFERENCE  AERODYNAMICS 

Again,  complete  generality  of  geometry  is  not  possible.  The  basic  fin  plan- 
form  is  assumed  to  be  a  trapezoid  with  root  and  tip  edges  parallel  to  the  free 
stream.  No  camber  is  considered.  Either  planar  or  cruciform  combinations  of 
fins  are  considered.  Two  sets  of  inline  lifting  surfaces  are  considered  in  the 
most  general  case.  Longitudinal  control  deflections  of  all  movable  surfaces  are 
assumed.  In  addition  to  the  individual  fin  geometry,  the  location  of  the  fin 
apex  and  a  mean  body  diameter  are  given  as  input  data. 

The  fin  cross  section  can  be  a  modified  double  wedge  or  biconvex.  Leading 
and  trailing  edges  may  be  blunted  cylindrically  and  independently  normal  to  the 
edges.  Each  piecewise  continuous  segment  of  a  fin  cross  section  s  assumed  to 
be  similar.  Figure  3  shows  a  typical  fin  planform. 
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Figure  3.  Fin  Planform  Geometry  for  a  Blunted 
Modified  Wedge  Cross  Section 


Planform  parameters  shown  which  are  needed  to  describe  the  modified  wedge 

are  y  ,  b,  C  ,  C.  ,  C  ,  and  C  .  For  a  blunted  biconvex  section,  C  and  C 

*1  ’  r’  t’  ri  r2  ’  rx  r2  are 

not  needed.  All  dimensions  are  in  feet.  The  x  and  y  are  local  coordinates. 

From  the  geometry, 

b/2  tan  y  +  C  -  C 

Un  \  = - V2 - -  •  (27) 

The  values  of  y  and  y  depend  on  the  cross-section  parameters.  The  dashed 
2  3 

lines  are  the  traces  where  the  cylindrical  bluntings  meet  inclined  planes  of  the 
wedge.  A  root  cross  section  for  a  modified  wedge  is  shown  in  Figure  4.  The 
slope  is  continuous  where  the  cylindrical  blunting  meets  the  inclined  plane. 


Figure  4.  Fin  Root  Cross  Section 
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Each  of  the  five  segments  shown  is  piecewise  similar  with  linear  taper  in  y.  In 
Figure  4,  the  radii  appear  as  segments  of  ellipses.  The  required  cross  section 
geometry  inputs  are  tf,  tt,  p^,  ptr,  p£t>  and  pu  where  tfc  is  the  tip  thickness, 
p^r  is  the  root  leading  edge  radius,  ptf  is  root  trailing  edge  thickness,  and  p^t 
and  ptt  are  similar  parameters  for  the  tip.  For  a  biconvex  cross  section  tt, 
p^t,  and  p  are  computed. 

From  the  planform  geometry, 

b/2  tan  Yi  + 

tan  Y  - - -  •  (28) 

2  b/2 


Here,  is  obtained  from  geometry  of  a  plane  perpendicular  to  the  leading  edge 
at  the  tip  and  the  assumption  of  piecewise  similarity. 

Cti  =  (tt/2  cos  6*  -  p£t  (1  -  sin  6i)]/(ain  6i  cos  Y^  (29) 

The  angle  is  the  angle  of  inclination  of  the  forward  plane  in  a  cross  section 
perpendicular  to  the  leading  edge  and  is  given  by 

sin  =  -  B  +VB2  +  C  (30) 


where 


p„_  (C  cos  Y  -  f  _) 


» '  K'W  *  nNrCK  fW-  ' 
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For  the  biconvex  case, 


plt  =  p£r  Ct/Cr  ’ 

*t  =  lr  Ct/Cr  1  and 
ptt  =  ptr  *VCr  ‘ 


(34) 

(35) 

(36) 


Equations  (34)  through  (36)  follow  from  the  requirement  of  piecewise  similarity. 


2. 1.2.1  AXIAL  FORCE 

Axial  force  components  are  broken  down  into  a  linear  superposition  of  wave 
friction  and  base  components  for  each  fin.  In  addition,  the  tail  affects  the 
incremental  base  axial  force  given  by  Equation  (20).  The  total  axial  force  for 
fins  is  given  by  the  sum  of  the  forces  per  fin  times  the  number  of  fins  (2  or  4 
assumed  per  set).  The  canard  and  tail  fin  set  forces  are  added.  The  only  inter¬ 
ference  is  the  incremental  base  force  due  to  the  tail  (four  fins  only).  In  addi¬ 
tion,  the  fin  normal  forces  in  the  presence  of  a  body  in  the  case  of  a  fin  de¬ 
flection  have  component  sin  6p.  contributions  to  the  axial  force.  The 

angle  6^  is  the  fin  control  deflection. 

2. 1.2. 1.1  WAVE  AXIAL  FORCE 


2. 1.2. 1.1.1  LOW  SUPERSONIC  FLOW  (1.05  <  M@  <  M^).  The  transonic  wave  axial 
force  is  a  simple  lirt.;r  descay  to  zero  from  the  value  at  -  1.05  given  by  the 
model  presented  here. 

The  basic  model  equation  is  the  first  order  linear  potential  model  which 
was  applied  to  >  he  first  order  lift  and  drag  problem  for  the  body-alone  case. 

P2  -  $  -  <ji  =0  (37) 

xx  Tyy  Tzz 

where 

<J>  is  the  disturbance  potential, 

x,  y  lie  in  the  plane  of  the  free  stream  velocity  vector,  and 
z  is  up. 
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The  boundary  conditions  are  applied  at  2  =  0  and  the  spanwise  contributions  are 
neglected.  A  source  distribution  in  the  x,y  plane  is  used  to  satisfy  the  dif¬ 
ferential  equation.  Application  of  the  limit  as  z  =  ±  0  determines  the  source 
distribution  strength  as  a  local  solution  in  analogy  with  slender  body  theory. 
Thu3,  for  the  lift  or  drag  problem, 


♦(x*y>t0)  ■ ;  "  f  10 


w(£».  n)  jL  lo_ 

f(K)2-P2(y-n)2 


where  a,  the  hyberoolic  influence  region  upstream,  is  defined  by 

(x-£)2  -  P2  (y-n)2  =  0  . 


(38) 


For  the  drag  problem, 

w(£>  n)  a  ±  |  f|  (€,  n) 


where  t  is  the  thickness  distribution.  From  similarity,  the  slope  distribution 

is  the  same  in  every  chord  plane  for  every  piecewise  continous  part  of  the  cross 

section.  For  the  biconvex  case,  (x/c)  is  continuous.  For  the  modified  wedge, 
a*.  GX 

is  discontinuous.  The  integration  region  is  bounded  by  the  projection  of 
the  upstream  Mach  cone  onto  the  x,y  plane  and  its  intersection  with  the  leading 
edge  and  side  edges  of  the  planform  consisting  of  two  finB  placed  together.  No 
intersection  of  the  opposite  side  is  considered.  The  numerical  method  of  Refer¬ 
ence  2  superimposes  basic  solutions  for  constant  slope  regions. 

The  0^  (x,y)  distribution  is  given  by 

cp  (x,y)  =  -  2  |J  .  (39) 

Note  that  this  is  the  pressure  coefficient  for  a  -  0,  =  0.  For  a  subsonic 

generator  (lines  associated  with  angles  y  ,  Y  ,Y  >Y  )>  C  has  an  integrable  singu- 

1  2  3  4  P 

larity. 

For  the  blunted  leading  edge  case,  a  modified  Newtonian  distribution  is 
assumed  to  apply  from  the  stagnation  point  to  the  end  of  the  cylindrical  section. 
The  pressure  coefficient,  C^,  is  then  discontinuous  at  this  "match"  point.  The 
thin  wing  theory  is  used  from  the  forward  junction  to  any  similar  junction  with  a 
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trailing  edge  radius.  It  is  assumed  that  separation  occurs  after  the  rear  cylin¬ 
der-airfoil  junction  and  that  no  contribution  to  wave  axial  force  occurs  past 
this  junction.  The  fin  base  force  is  considered  later. 


Once 

given  by 


L  is  obtained,  the  wave  axial  force  coefficient 
P 


b/2  ytan 


Y  +  C(y) 

l 

(x,y)  dxdy 


for 


0  ytan  Y 

l 

Here  C(y)  is  the  local  chord.  For  the  wing  alone  case  S  ^  is 
for  two  fins.  Otherwise, 


a  single  fin  is 


(40) 


the  planform  area 


ref 

Note  that,  for 
Y  =1°.  For 

4 

given  planform. 


Y  <  0  or  y  <  0,  the  planform  is  modified  so  that  Y  =  1°  or 
3  4  3 

such  a  planform,  the  resultant  axial  force  is  higher  than  for  the 

For  greater  detail,  see  Reference  2. 


2. 1.2. 1.1. 2  TRANSONIC  FLOW  (M^  <  1.05).  A  linear  decay  from  the  value  at  = 

1.05  to  zero  at  If  -  .85  is  assumed  such  that 
00 


CAW  =  5  CAW(M  =1.05)  (M« 
00 


.85),  .85  <  <  1.05  . 


(41) 


CAW  =  0  for  <  -85  • 


2. 1.2. 1.1. 3  HIGH  SUPERSONIC  FLOW  (M^  >  M^).  The  theory  used  to  determine  the 
lifting  surface  properties  is  given  and  evaluated  in  Reference  4.  The  theory 
used  to  determine  the  pressure  is  based  on  two-dimensional  supersonic  flow  prop¬ 
erties.  At  a  given  point  on  the  surface,  the  slope  of  the  surface  relative  to 
the  undisturbed  free  stream  i"  determined  first.  If  this  angle  makes  a  compres¬ 
sion  surface,  the  pressure  is  calculated  from  oblique  shock  theory  using  the  in¬ 
clination  angle.  If  the  angle  indicates  an  expansion  surface,  the  pressure  is 
calculated  from  Prandtl-Meyer  theory.  There  is  no  spanwise  variation  of  the 
pressure  distribution.  For  blunt  leading  or  trailing  edges,  a  Newtonian  distribu¬ 
tion  is  assumed  as  before. 
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The  axial  force  is  given  by  an  equation  similar  to  Equation  (40).  How¬ 
ever,  in  this  case  the  lift  and  drag  problems  are  solved  simultaneously.  The 
single  fin  axial  force  is  given  by 


b/2  C(y)  +  ytan  y 


C  =  — ■ 
2S 


J  / 


CpU  +  CpL/3x  dxdy 


ytan  y 


where  U  and  L  refer  to  upper  and  lower  surface  values,  respectively. 

2. 1.2. 1.2  SKIN  FRICTION  AXIAL  FORCE 

The  model  here  is  essentially  the  same  as  for  the  body  alone  case.  For  a 
single  fin, 


CAF  = 


where 


-  »  /l. 328  .  - 

F  Wc  K 


S  »  b(C  +  C. )/2  and  is  the  wetted  area  for  a  single  fin. 
w  r  t 

Here,  R^c  is  5  x  10s  and  is  based  on  the  mean  aerodynamic  chord 

(c?  -0 

n  —  o  /o  > _ r—  A  L 


C  =  2/3 


(C?  '  Cr  ) 


From  flat  plate  subsonic  strip  theory,  the  fin  x^  lies  on  the  quarter  chord 
where  C  =  C. 


2. 1.2. 1.3  TRAILING  EDGE  SEPARATION  FORCE 

The  model  is  based  on  the  base  pressure  for  flow  over  a  two-dimensional  or 
semi-infinite  (in  spanwise  direction)  flat  plate.  This  pressure  is  given  by  a 
table  lookup  interpolation  function  of  Mach  number.  For  a  single  fin, 

b  WH»> 

CAB  2S  'ptr  +  pttJ  ' 
ref 
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fHP 


2. 1.2. 1.4  EFFECT  OF  A  FIN  ON  THE  BASE  DRAG  OF  THE  BODY 

In  this  case,  Equation  (20)  is  replaced  by 

ACABA  =  1“  (-0035  "  -01  XAFT^  +  CpB 

-  .lx^)]  (RB/Href)3  .  (46) 

The  last  part  of  Equation  (46)  is  dropped  for  t^/C^  <  •  lx^pp  wliere  xaft  is  the 
distance  from  the  base  to  the  root  trailing  edge  divided  by  the  root  chord.  For 
flares,  (Rg/R^^)2  is  used  in  Equation  (46). 

2. 1.2. 2  LIFT  AND  PITCHING  MOMENT 

Lift  contributions  are  broken  down  by  the  component  superposition  principal. 
Here,  for  the  wing  alone,  which  consists  of  two  fins  placed  together,  wing  load¬ 
ing  is  considered  first.  Subsequently,  interference  effects  are  considered. 


2. 1.2. 2.1  SUBSONIC  FLOW  (0  <  M^  <  .8) 

Equation  (37)  applies  again.  The  well-known  Gothert  rule  transforms  the 
equation  to  Laplaces  equation  for  subsonic  Mach  numbers.  From  the  well-known 
solution  for  two-dimensional  incompressible  flow  over  a  cylinder,  a  bound  cir¬ 
culation  may  be  added  to  the  solution.  The  circulation  contributes  the  lift.  An 
integral  representation  of  the  three-dimensional  vector  may  be  obtained  by  de¬ 
fining  a  vector  potential.  The  vector  potential  satisfies  a  vector  Poisson  equa¬ 
tion.  This  leads  to  the  well-known  analogy  with  the  Biot-Savart  law  for  electro¬ 
magnetic  theory.  A  Kutta  condition  is  applied  at  the  trailing  edge  to  remove  the 
indeterminancy.  After  much  manipulation  and  application  of  boundary  conditions, 
it  may  be  shown  (References  19  and  20)  that 


p  tan  a  «  p  a 


8n 


f 


b/2 


-b/2 


^  XTE0(y) 

xLE0(y) 


ACp0  (X0’Y) 

(Fy)*' 


1  + 


(xQ  -  X0)2  +  (y 


Y): 


dXQdY 


(47) 
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where 


A  CpO  =  CpLO  "  CpUO  =  2V0.  P  =VT-  Kj  (48) 

and 

^  CpO  t^le  differential  pressure  coefficient  loading,  and 

Y0  is  the  streamwise  circulation  per  unit  length. 

The  outer  integral  requires  a  Mangier  principal-value  treatment  (Reference  20). 
Here,  x0  and  XQ  are  the  stretched  dimensions  x/p,  X/p.  The  pressure  coefficient 
differential  is  given  by 

ACp0  * 

where  the  subscript  0  refers  to  the  analogous  M„,  =  0  case.  The  compressibility 
parameter  p  cancels  when  AC^  is  substituted  into  Equation  (47). 

Equation  (47)  is  first  transformed  to  dimensionless  variables 


*0  •  V*)  xo  -  *le0«) 

X  *TP  (y)  -  “re  (y)  ’  £  =  x-„  (Y)  -  x. 


rr  ,  and 


t  =  JL  n  .  JL 

^  b/2  ’  "  b/2  * 

Substituting  Equations  (49)  and  (50)  into  Equation  (47)  yields 

a(x’y)  =  m  fo*c  p  n)  K(£*  n)d^ dn 


(49) 

(50) 

(51) 


K(C,n)  =  1  + 


(xo  _  V2  +  (y  '  Y>2 


(53) 


Equation  (53)  is  discontinuous  in  XQ  at  xQ  when  Y  =  y.  When  XQ  <  xQ,  K  =  2. 
When  XQ  >  xQ,  K  =  0.  a(x,y)  =  a  for  the  static  lifting  problem.  AC  is  given 
by  an  approximating  function  with  undetermined  coefficients.  Equation  (51)  is 
evaluated  once  for  each  of  the  coefficients  (method  of  collocation).  The  coef¬ 
ficients  are  obtained  by  matrix  inversion. 
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An  assumed  functional  form  (see  Reference  21)  is  given  as 


where 


n  m 

:  (£.  n)  =  E  £  Ac  a  ,  nJh  (fjg  (n) 

p  j=l  i=l  p  1  J  1  J 


(54) 


(55) 


8j(n) 


.  Hl oL 
V’-lnjI 


n  (|n|  -ln„|) 

m _ L- 

n  (|n,|  -  I  nj  ) 


The  interpolation  pivot  points,  (£^,  4^),  are  given  by 


h  =  i  1  ■ cos 


Values  of  Xri  chordwise  location  where  the  boundary  conditions  are  applied,  are 


Xr  =  1  "  4i»  r  =  !»  2.  . . .  ro,  i  =  m  - 


r  +  1  . 


The  interpolation  points,  r\.f  are  given  by  the  n  zeros  of  the  Jacobi  polynomial 

J 

n  (2m+2J2+l) !  qf 

M(nj)  =  °  =  jJb  2^  W  (m-£)f  (m+£) !  (2m+l)!  *  (59) 


The  boundary  condition  points  £  ,  s  =  l,  2,...  n  are  the  same  as  the  r|..  The 

®  J 

boundary  conditions  are  thus  evaluated  at  a  set  of  n  x  m  punts  (x  ,  £  )  r  =  1, 

r  s 

2,...  m;  s  =  1,  2,...,n.  The  functional  form  for  AC^  contains  the  leading  edge 

singularity  for  a  two-dimensional  flat  plate,  satisfies  the  Kutta  condition 

(ACp  =  0)  at  the  trailing  edge,  and  satisfies  the  tip  loading  condition  (AC^  =  0) . 

The  program  is  coded  effectively  as  m  =  3  and  n  =  8.  However,  four  r|  points 

are  chosen  for  r|  >  0  and  four  points  are  chosen  for  r)  <  0.  Two  functions  aie 

thus  needed  for  AC^,  one  for  r|  >  0  and  one  for  r|  <  0.  For  the  lifting  problem, 

AC  (£.,  n.)  is  even  valued  in  n. 

P  *i’  'j 
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The  normal  force  and  pitching  moment  (xCg  =  0)  coefficients  are  given  by 


b/2  ytan  y  +  C(y) 


C  s  — 
N  S 


/  /  v- 


y)  dxdy  ,  and 


ytan  y 


S  ,  C  , 
ref  ref 


b/2  ytan  y^  +  C(y) 
/  /  <x  +  “  Cre£> 


ACpdxdy 


ytan  y 


where  x  is  the  location  of  the  wing  apex  from  the  origin  in  calibers.  For  the 
wing  alone  case,  x  =  0  and  =  C  as  given  in  Equation  (44),  and  Srg^  is  the 
planform  area  for  two  fins. 


2. 1.2. 2. 2  LOW  SUPERSONIC  FLOW  (1.2  <  <  M£  ) 

e 

The  Mach  number,  M^  ,  separates  low  supersonic  computations  from  high  super¬ 
sonic  computations  for  tRe  wing  static  and  dynamic  lifting  and  interference  aero¬ 
dynamics.  M^  is  a  program  input  chosen  as  high  as  practical  for  the  thin  wing 
theory  application  (usually  4.0  <  Mg  <  5.0)  as  this  provides  a  better  comparison 
with  inviscid  lift  and  moment  coefficient  data  than  the  strip  theory. 

The  basic  solution  of  Equation  (37)  for  a  thin  wing  is  given  by 

♦  «;i  f  f  n)di  an  ,62) 

"  Jo  JT|c*-4>2-|52(y-n)2 

where 

w(£,  r))  =  *  a  on  the  wing  only. 

In  this  case,  the  area  of  influence,  o,  contains  points  which  are  not  on  the 
wing.  For  points  off  the  wing,  w(£,  r|)  is  not  zero  as  was  the  case  for  the  sym¬ 
metric  axial  force  problem. 

The  leading  edge  may  be  supersonic  or  subsonic.  The  trailing  edge  is 
restricted  to  the  supersonic  case.  If  the  trailing  edge  is  subsonic,  the 
input  Mach  number  is  increased  so  that  the  trailing  edge  is  slightly  super¬ 
sonic.  The  swept-forward  trailing  edge  case  presents  no  problems.  Again  no 
influence  of  an  opposite  half  wing  tip  is  assumed;  the  Mach  number  is  increased 
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for  those  low-aspect  ratio  planforms  where  this  occurs.  For  a  more  detailed 
explanation  and  evaluation,  see  Reference  2. 

The  supersonic  leading  edge  case  is  the  simplest.  Three  subcases  are  de¬ 
lineated  . 

Figure  5  indicates  Case  1  when  the  Mach  lines  from  the  apex  and  leading  edge 
at  the  tip  do  not  intersect.  For  this  case,  three  zones  bounded  by  the  x  axis, 
right  Mach  line  from  the  root  apex,  and  left  Mach  line  from  the  tip  apex  are 
defined.  The  regions  are  designated  as  1,  2,  and  3.  Typical  a  areas  of  in¬ 
tegration  are  bounded  by  correspe tding  numerals  1,  2,  and  3.  For  region  3,  the  a 
area  of  influence  is  determined  by  the  method  of  Eward  and  Krasilchikova  (Refer¬ 
ence  20).  For  explicit  ACp  relations,  see  Reference  2. 


Figure  5.  Supersonic  Leading  Edge  Flat  Plate  Planform,  Case  1 

<M„  =  2) 

A  second  supersonic  leading  edge  situation,  Case  2,  is  shown  in  Figure  6. 
Here  the  two  Mach  lines  intersect  and  the  root  chord  apex  Mach  line  intersects 
the  trailing  edge.  Mw  is  lower  fee  this  case. 
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9  - - -  DIVIDING  MA*1M  LIN££  POR 


Figure  6.  Supersonic  Leading  Edge  Flat  Plate  Planform, 

Case  2  (M^  =  1.5) 


Four  zor,?s  are  obtained  this  time,  as  indicated  in  Figure  6. 

In  Figure  7  (Case  3),  the  apex  Mach  line  intersects  the  side  edge.  Only 
the  integration  area  for  zone  5  is  shown  this  time.  In  this  case, 

'.  ACp  =  (ACp)2  -  (ACp)s. 

(ACp)c,  refers  to  an  equivalent  integration  tip  region  contribution.  The  approach 
taken  in  Reference  2  seems  to  be  more  complex. 

The  approach  for  the  subsonic  leading  edge  case  is  much  more  complex  than 
the  supersonic  leading  edge  case.  Figure  8  indicates  the  regions  of  interest  and 
integration  influence  region. 

For  the  case  of  Figure  8,  no  part  of  the  side  or  leading  edge  is  supersonic. 
Here,  w  is  unknown  between  the  planform  and  the  apex  Mach  lines.  For  region  1, 
the  solution  is  given  by  transforming  Equation  (37)  in  conical  flow  variables  to 
a  two-dimensional  Laplace  equation  in  a  cross  section  perpendicular  to  the  Mach 


cane  from  the  apex,  for  points  in  region  2,  an  additional  tip  solution  must  be 
added  to  cancel  lift  off  the  wing.  The  AC^  relations  are  given  in  Reference  2. 

2. 1.2. 2. 3  TRANSONIC  FLOW  (.8  <  M  <1.2) 

00 

Thickness  effects  are  neglected  for  the  subsonic  and  supersonic  flow 
regimes.  The  approach  taken  in  transonic  flow  is  the  same  as  that  in  DATCOM 
(Reference  22).  Details  are  also  given  in  Reference  2. 

A  linear  interpolation  in  a  table  of  Cjj  values  is  used.  At  =  .8  and 
M  ~  1.2,  Cu  is  computed  using  the  subsonic  and  low  supersonic  models.  Three 

"ot 

additional  MJJch  numbers  are  then  computed:  Mfb,  Mfl  =  M^b  +  .07,  and  M^  = 

+  .14  where  Mfb,  the  force  break  Mach  number,  is  obtained  by  a  table  lookup  in¬ 
terpolation.  First.,  a  force  break  Mach  number,  (M-^)o*  for  an  equivalent  rec¬ 
tangular  planform  is  obtained  as  a  table  lookup  function  of  t/C  and  aspect  ratio. 
The  thickness  t  is  taken  at  the  mean  aerodynamic  chord  C  given  by  Equation  (44) . 
Then.  M^  is  obtained  as  a  table  lookup  function  of  (M^)q  and  the  sweep  of  the 
mid-chord  line,  Next>  CN  at  Ma»  and  ^  are  computed  where  the  the¬ 

oretical  value  of  CN  at  Mfb  is  ¥irst  given  from  lifting  line  theory  as 


(V£bt 


2n(AR> 


2  +  [ARZ(P2  +  tan2  Y1/2)  +  *1 


AR  is  the  aspect  ratio,  and 


f>  4h  ■ 


Then , 


(cu  ) 

cr  )  =  (C  )  »  £b 

V£b  V£bt  (CK  >fbt 

a 


Tha  r.'itio  in  Equation  (64)  is  given  as  a  table  lookup  function  of  t/C  and  aspect 
ratio.  Then  (CM  ■.  and  (CM  ),  are  given  as 

Va  a  D 

(CN  \  =  (1  '  Ka}  (CN  }fb  ’  and  (65) 

a  a 


(66) 


(CN A  =  (1  “  V  (CN  Jfb 

a  a 

where  is  given  as  a  table  lookup  function  of  t/C  and  AR  and  is  given  as  a 
table  lookup  function  of  t/C.  The  final  value  of  Cjj  is  obtained  by  a  linear 
interpolation  in  a  table  of  the  five  values  which  have  Been  computed. 

The  x£  is  given  by  a  linear  interpolation  between  the  value  at  Ma  s  .8  and 
the  value  at  =  1.2. 


2. 1.2. 2. 4  HIGH  SUPERSONIC  FLOW  (Mtf  >  ) 

e 

The  tangent  wedge  method  has  been  previously  outlined  in  connection  with  the 
wing  wave  axial  force, 
one  obtains: 


For  lifting  properties,  defined  in  terms  of  and 


b/2  ytan  y  +C(y) 


CN  =  S 


ref 


/  /  V 


C^)  dxdy  ,  and 


(67) 


ytan  y 


b/2  ytan  y  +C(y) 


r  r  -  _ - _ 

m  S  ,  C  , 
ref  ref 


If  “:i 

0  ytan  yi 


Cref)(CpL  '  V  dxdy  * 


(68) 


2. 1.2. 3  INTERFERENCE  LIFT 

The  total  lift  as  given  by  component  buildup  is  discussed  rather  extensively 
in  Reference  2.  For  ,  interference  is  not  considered.  For  ,  in¬ 
terference  factors  are  obtained  from  slender  body  with  certain  correction  factors 
taking  into  account  canard  and  body  vortex  wakes  as  well  as  lifting  surface 
proximity  to  the  body  base.  Compressibility  corrections  are  also  made. 

2. 1.2. 3.1  LIFTING  SURFACE-BODY  INTERFERENCE 

The  first  problem  to  be  considered  is  that  of  an  undeflected  lifting  sur¬ 
face-body  combination  at  aagle-of-attack.  For  the  combination,  the  lift  is  given 
by 


.32 


where 


Lc  =  Kc  lf 


(69) 


Kc  =  KB(F)  +  KF(B)  +  kb  * 

(70) 

KB(F)  =  LB(F)  1  **  ’ 

(71) 

^(B)  =  LF(B)  1  LF  *  and 

(72) 

kb  =  lb  '  h 

(73) 

and  where 

Lp  is  the  lift  of  the  fin  alone  (2  fins  placed  together), 

Lb(F)  on  t*le  due  to  the  presence  of  the  fins, 

Lp(B)  is  the  lift  on  the  fin  in  the  presence  of  a  body, 

Lg  is  the  body  alone  lift, 

Lg  and  are  predicted  from  the  theories  already  described,  and 

and  are  given  by  slender  body  theory  (see  Reference  23). 

Figure  9  shows  a  triangular  fin  and  body  combination  planform  view. 


Figure  9.  Triangular  Fin  and  Body  Combination 


Slender  fin-body  theory  (radial  distance  small  from  the  axis)  shows  that  the 
lift  problem  solution  is  given  by  a  local  solution  of  Laplaces  equation  in  the 
cross  flow  plane.  Solution  of  Laplaces  equation  is  readily  obtained  with  the 


aid  of  the  Joukowski  and  Swartz-Christoffel  transformations.  R  is  a  mean  local 
radius . 


The  resulting  vertical  loading  coefficients  are  given  in  terms  of  R, 

Y,  and  y.  Integration  over  the  individual  wing  and  body  planform  areas  gives 

If  -  4  +  ^2)  „/9  (a  +  a)  -1/  2A  \ 

^(B)  "  tt(A-TP  [  £  2  4  Ur5rl/ 


^  -  j)1, 

2 


W(A) ,  and 


KB(F)  =  (!  +  a)  ‘  KF(B)  = 


(74) 

(75) 


where 


A  =  b/2  +  R 
R 


If  the  trailing  edge  remains  straight  or  swept  forward,  the  results  in 
Equations  (74)  and  (75)  apply.  However,  for  the  swept  forward  trailing  edges, 
KF(B)  ancl  Kjj(F)  are  zero  if  the  pressure  loading  is  integrated  to  the  trailing 
edge.  The  argument  given  in  Reference  23  is  that  the  planform  beyond  where  y 
is  a  maximum  contributes  nothing  to  the  loading  since  the  forward  downwash 
cancels  the  angle-of-attack  loading  on  the  rear.  Therefore,  the  triangular 
planform  results  are  assumed  to  hold  for  the  trapezoidal  planform. 

For  the  swept  back  trailing  edge  in  Figure  10,  the  solid  boundaries  for 
the  cross  section  are  different  downstream  from  the  root  trailing  edge.  In 


Y 


Figure  10.  Swept  Back  Trailing  Edge  Fin  and  Body  Combination 
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Reference  2  the  following  modifications  are  made: 


KB(F)  = 

B  (A)  G  ,  and 

(76) 

KF(B)  = 

1  +  (W  (A)  -  1)G 

(77) 

where 

G  ‘  b/2 

C 

r 

tan  y  + 

(78) 

When  B(A) 

=  0,  W(A)  =  1,  the  large  aspect  ratio  limits  are  obtained. 

The  G  =  1 

limit  corresponds  to  straight  or  swept  forward  trailing  edges.  The  case  when 
b/2  tan  <  cr  <  b/2  tan  Yj  +  Ct  needs  additional  investigation.  For  this  case, 
the  cross  sections  are  the  same  up  to  the  root  trailing  edge  as  that  for  which 
Equations  (74)  and  (75)  were  derived.  The  forward  downwash  would  then  cancel 
the  rear  planform  loading.  For  the  case  when  <  b/2  tan  Yj »  G  might  be  re¬ 
placed  by 

(b.)/2  +  R 

G  =  - -  . 

b/2  +  R 

In  addition,  thete  is  a  lift  contribution  by  a  finite  control  deflection 

CNF(B)  =  ^(B)  ^N^F  6F  ‘  (79) 

The  potential  problem  in  the  crossflow  plane  (a  =  0)  consists  of  a  uniform  down- 
wash  on  the  lifting  surface  in  an  undisturbed  far-flow  field  (see  Reference  23). 
In  addition, 

CNB(F)  =  (KF(B)  "  ^F(B)^  (CN  >F  6F 


=  (W(A)  -  A(A))  (CN  )F  6f  .  (80) 

An  explicit  relationship  for  is  given  in  Reference  2.  For  a  swept  trail¬ 

ing  edge  a  correction  is  made  as  before: 

kF(B)  =  1  +  (A  (A)  -  1]  G  ,  and  (81) 

kB(F)  =  [W(A)  ‘  A(A)]  G  ‘  (82) 

Equation  (75)  for  B  is  used  for  subsonic  Mach  numbers  and  Mach  numbers  up  to 
when  the  Mach  line  from  the  leading  edge  of  the  tip  does  not  intersect  the  body 
section  perpendicular  to  the  root  chord.  Figure  11  illustrates  three  possible 
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:igure  11.  Determination  of  Kg^  for  High 'Aspect 
Range  at  Supersonic  Speeds 


subcases.  In  this  case,  the  body  loading  area  is  redefined  by  the  crosshatching. 
The  criteria  for  this  case  is  given  by 


2C 

P  >  "  tan  *1 


(83) 


For  the  no-afterbody  case,  is  a  function  of  C^/0.  For  0D  >  (the  Mach 
line  from  the  root  chord  leading  edge  intersects  the  end  of  the  body),  C^/0 
is  set  equal  to  D.  Subcases  for  the  no  a^ter  body  case  are  subsonic  and  super¬ 
sonic  leading  edges.  Subsonic  leading  edge  values  of  are  given  for  the 
infinite  afterbody  value  case  of  x^pp.  The  short  afterbody  computation  of  Bg^ 
is  given  by  the  linear  interpolation  where  (x^pp  <  D0) : 
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2. 1.2. 3. 2  WING-TAIL  INTERFERENCE 


Thus  far  only  body-lifting  surface  interferences  have  been  considered.  An 
additional  interference  is  given  by  vortex  effects  of  the  canard  or  wing  on  the 
tail.  A  single-line  vortex  is  assumed  to  be  shed  from  each  canard  panel.  The 
vortex  strength  relationship  and  effective  lateral  location  are  obtained  from  the 
wing  loading  distribution  given  by  slender  wing-body  theory.  The  vortices  are 
assumed  to  be  aligned  with  the  free  stream  velocity.  The  canard  hinge  is  assumed 
to  be  at  the  quarter  chord.  This  is  enough  information  to  adequately  define  the 
vortex  locations  at  the  tail.  The  negative  lift  coefficient  due  to  the  downwash 
of  the  canard  shed  vortex  on  the  tail  is  given  by 


CNT(V)  =  (CN  \  (CN  ^(B)  Sin  “ 
a  a 


+  “W(B)  Si“  V 


ip  (bp  -  Rp  ) 


3W 


2n  ARp  (fw  -  Ry)  \  ref 


(85) 


Here,  the  W  subscript  refers  to  the  wing  or  canard,  the  T  subscript  refers 
to  the  tail,  f^  is  the  span  location  of  the  panel  vortex  obtained  from  slender 
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body  theory,  and  i,p  is  the  tail  interference  factor  also  obtained  from  slender 
body  theory.  The  crossflow  potential  problem  for  the  tail-body  combination  con¬ 
sists  of  the  solution  for  two  point  vortices  of  a  given  strength  and  location  in 
an  undisturbed  flow  field.  b(^  is  the  tail  half-span  (including  the  body  radius). 

In  Reference  2  an  expression  is  given  for  a  negative  lift  on  the  tail 
due  to  body  shed  vorticity.  This  is  neglected  in  the  current  version  of  the 
program. 


The  normal  force  coefficient  can  now  be  written  as 


CN  ~  CNB  + 


*<[= 


{[*» 


(B)  +  KB(W)]  “  *  [kW(B)  +  kB(W) 


KT(B)  +  KB(T)]  a  +  [kT(B)  +  kB(T 


]  M  VW 

)]6T  }(CNa}T 


+  C 
=  CNB 


NT(V) 
+  C 


NW(B)  +  CNB(W)  +  CNT(B)  +  CNB(T)  +  CNT(V)  ’ 


(86) 

,  and 


At  subsonic  Mach  numbers  (M^  <  .8),  the  center  of  pressure  for 
CNB(W)  are  assume^  to  be  at  the  win8  alone  center  of  pressure.  For 
CNT(B)’  and  CNB(T)’  the  Xcp  are  assuaiecl  to  coincide  with  the  tail  alone  x  . 


For  supersonic  flow  CNfi^  and  cjjb(x)  * 


the  x  values  are  assumed  to  be 
cp 

given  by  the  centroid  of  the  areas  in  Figure  11.  At  transonic  speeds  the  center 
of  pressure  values  are  given  by  linear  interpolation  between  the  values  at 
M  =  1.2  and  M  =  M, 


00 


00 


fb" 


2.1.3  HIGH  ANGLE-OF-ATTACK  AERODYNAMICS 


The  current  approach  is  empirically  based  on,  and  is  an  adaption  of,  the 
method  described  in  Reference  24. 


The  range  of  input  parameters  for  a  cruciform  configuration  is  limited  to: 

1.  .8  <  M  <  3.0; 

00 

2.  Angle-of-attack  0°  to  180°  for  isolated  components  (roll  angle  0°)  and 
0°  to  45°  for  body-tail  combinations  at  arbitrary  roll  angles  from  0°to  180°; 


38 


3.  Tail:  trapezoidal  planform,  side  edges  parallel  to  the  body  centerline 

A.  Leading  edge  sweep  angle:  0°  to  70°; 

B.  Taper  ratio:  0  to  1; 

C.  Aspect  ratio  (two  panels):  .5  to  2.0; 

D.  Trailing  edge  straight  and  parallel  to  end  of  body; 

4.  Nose  length  (pointed  tangent  ogive):  1.5  to  3.5  calibers; 

5.  Cylindrical  afterbody:  6  to  18  calibers; 

6.  Total  span-to-diameter  ratio  (two  fins):  1  to  3.33. 

The  roll  definition  and  positive  fin  load  orientation  are  shown  in  Fig 
ure  12.  Note  that  6  is  measured  from  the  windward  side. 


Figure  12.  Roll  Angle  and  Fin  Load  Definitions  (Looking  Forward) 


The  axial  force  coefficient  is  assumed  to  be  contributed  entirely  by  the  body. 
The  total  normal  force  coefficient  is  given  by 

CN  =  CNB  +  t(CNTl  +  SiTS^  Sin  9  +  (CNT2  +  ST45  C0S  9^  ST/Sref 
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BT 


(87) 
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where  are  the  individual  fin  loads  in  the  presence  of  the  body  and  other 

fins,  and  Ig^,  is  the  tail-to-body  carryover  normal  force  coefficient  ( i.e 

^NB(T)^  * 

The  longitudinal  center  of  pressure  from  the  nose  for  the  entire  configura¬ 
tion  is  given  by 

xcp  =  {xcpB  CNB  +  ^Sln  6  (xcpTl  CNT1  +  XcpT3  CNT3) 

+  cos  0  (xcpT2  Cm  +  xcpT4  ^j,T4)  ]  Sx/Sref  +  JBT  Xcpi  } /CN  (88) 

where 

x  values  are  given  in  calibers, 

x  D  is  the  body  alone  x  , 
cpB  7  cp’ 

XcpTi  are  the  lndividual  tail  xcp  values  iu  the  presence  of  the  body  and 
other  fins,  and 

xcpj  is  the  center  of  pressure  for  the  carryover  load. 

if  one  neglects  the  small  body  frictional  roll  coefficient  and  the  tail 
carryover  onto  the  body,  one  can  estimate  the  static  roll  coefficient  as 

CJ£  =  2  bVSref  X.  CNTi  S1U  ”  1  *  ^  ^ 

i=l 

+  b/D  iyrpi  ;  b/2)]  •  C89) 

Details  of  the  functional  forms  are  given  in  Reference  24. 

2.2  DYNAMIC  DERIVATIVES 

Dynamic  derivative  prediction  is  for  email  angles-of-attack  (a  0)  only. 

2.2.1  BODY  ALONE 

2. 2. 1.1  ROLL  DAMPING  AND  BODY  MAGNUS  MOMENT 

These  coefficients  are  given  empirically.  The  data  base  is  configuration 
limited.  Coefficient  Mach  number  dependence  was  0  <  <  5.  An  asymptote  is 
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assumed  at  M^  =  5.  Geometric  parameters  entering  the  empirical  functional  form 
are  the  total  length,  center  of  moments  location,  and  boattail  length,  all  in 
calibers.  The  boattail  length  is  limited  to  1.5  calibers  for  the  Magnus  moment 
coefficient  computation.  The  Magnus  moment  coefficient  C  ^  is  computed  at  1° 
and  5°  angles-of-attack.  Magnus  moment  coefficients  are  not  estimated  for  bodies 
with  lifting  surfaces. 
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2. 2. 1.2  PITCH  DAMPING  COEFFICIENT 

The  original  prediction  used  the  empirical  prediction.  The  method  described 

in  References  25  and  26  was  intended  to  improve  the  C  +  C  •  prediction.  Evaiu- 

mq  ma 

ation  of  the  method  is  given  in  Reference  4.  Restrictions  are: 

(1)  Initial  cylindrical  radius  for  a  spherically  blunted  or  trun¬ 
cated  body  is  limited  to  less  than  .25  calibers. 

(2)  Contribution  of  afterbody  and  boattail  or  flare  are  neglected 
at  higher  Mach  numbers. 


In  subsonic  flow,  C  +  C  •  is  given  by  a  relation  based  on  slender  body 

mq  ma  °  1 


theory: 


c  +  C  •  =  -  4  (-77  +  .23  M^)  CL  V-  77  +  .23  M?  -  x  )2  , 
mu  ma  <»  y  00  ex 

■& 


(90) 


mq  ma  »  ’  00  eg 

Here,  L  is  the  total  length  in  calibers,  and  X(<^  is  the  center  of  moment  location 
from  the  nose  in  calibers. 


Hypersonic  (embedded  Newtonian)  flow  approximations  are  applied  .  •o\'e  a 
certain  Mach  number  M  =  M*  associated  with  an  effective  hypersonic  >  itity 
parameter  of  .4  such  that 


M*  =  .4  esc  0* 
-1  /  1 


0*  =  tan 


JN 


(91) 

(92) 
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If  M*  <  1.5,  then  M*  =  1.5  is  chosen.  For  1  <  II  <  M*  a  linear  inter- 
7  00 

polation  is  used  where 


M*  - 

Snq  +  ^raa  ~  M*  -  1  ^mq  +  Mw  =  1 

M»  "  1 

+  M*  -  I  ^mq  +  ^m<P  Mw  =  M* 


(93) 


For  M^  >  M*,  a  direct  application  of  a  modified  Newtonian  pressure  dis 
.ribution  at  a  =  0  for  the  pitching  motion  gives 


C  +  C  *  =  -  16 
mq  ma 


/» 

0 


C  C  (x  -  x  +  RR')2  RR’ 

-£9— I.- _ _£S_ _ _ _  dx 

1  +  R  1  * 


(94) 


Note  that  for  a  spherically  blunted  body  the  contribution  of  the  spherical  cap 
is  neglected.  The  stagnation  pressure  coefficient,  C^o,  behind  a  bow  shock  is 
approximated  for  large  Mach  numbers  as 


=  ULA 
Y  +  i 


1.5 
Y  +  3 


(95) 


and  C  is  given  as 

3 

I  1.01  +  1.31  [log  (10  sin  e)]"7/3  ]  sin  e  >  .4 
C  =  (96) 

'  I  1.625  ) M  sin  e  <  .4  . 

The  Lockheed  Missiles  and  Space  Company  (LMSC)  model  was  assumed  to  be  able 
to  replace  the  empirical  method  for  the  body-alone  pitch  damping  coefficient. 
Since  the  LMSC  prediction  neglects  the  effects  of  afterbody  and  boattail  at 
supersonic  Mach  numbers,  the  empirical  method  often  predicts  better  values  for 
pitch  damping.  For  this  reason  and  because  of  other  limitations  for  bodies  with 
lifting  surfaces,  the  pitch  damping  prediction  method  is  currently  an  input 
option  to  the  code. 

The  LMSC  body-alone  algorithm  was  extensively  evaluated  by  R.  Whyte  in  Ref¬ 
erence  27.  A  conclusion  reached  in  Reference  27  is  that  the  LMSC  algorithm  is 
adequate  for  most  applications.  For  certain  nose-afterbody  configurations  with 
the  located  forward,  poor  results  were  obtained  when  compared  with  data  at 


42 


supersonic  Mach  numbers.  The  code  was  modified  for  these  cases.  When  the  value 


of  C  predicted  is  less  negative  than  that  predicted  by  the  empirical  method  and 
not  within  75  percent  of  the  empirical  value,  the  empirical  value  is  chosei. 

2.2.2  LIFTING  SURFACE  DYNAMIC  DERIVATIVES 

The  rolling  rate  is  assumed  small  and,  hence,  Magnus  moment  is  neglected. 

2. 2. 2.1  ROLL  DAMPING 


2. 2. 2. 1.1  SUBSONIC  FLOW  (0  <  M  <  .8) 

Equation  51  is  again  applicable.  The  effective  angle-of-attack  is  given  by 

a(x,y)  =  ^ .  In  this  case,  the  resultant  ('AC  )  obtained  by  matrix  inversion  is 
00  I*  P 

odd-valued  in  r). 

With  CACp)^  =  (ACT)p  P/V*»  C£p  is  «iven  as 

%  =  C1  '  "2-^  (97) 

where  Crpf.p  is  b  for  the  wing-alone  case.  Using  the  definition  in  Equation  (97) 
one  obtains 


b' 
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s  _  c2  .  J  . 

ref  refp  0  ytan  y 


ytan  y  +  C(y) 
y(ACT)p  dxdy  . 


Note  that  the  planform  considered  is  that  obtained  by  extending  the  leading  and 
trailing  edges  to  the  body  centerline,  b'  is  the  total  span. 


2. 2. 2. 1.2  LOW  SUPERSONIC  FLOW  (1.2  <  M^  <  M£  ) 

e 

The  approach  used  for  the  static  lift  problem  is  applicable  here.  In  this 

case,  w  (4,  rj)  =  “  is  not  a  constant. 

00 

The  approach  for  subsonic  and  supersonic  leading  edges  is  similar  to  that 
used  for  the  static  ]ift  problem.  Specific  relationships  for  (ACp)p  _  p/V^  (ACp)p 


A3 


are  given  in  Reference  3.  The  modified  pianform  considered  for  subsonic  flow  is 
used  again. 


2. 2. 2. 1.3  TRANSONIC  FLOW  (.3  <  Jfl  <  1.2) 

Ol* 

An  empirical  approach  is  taken  here.  It  is  assumed  that  the  variation 
is  similar  to  that  for  C.T  such  that 

“  (CN  3M 

C  V «»  =  '  V«r.f  tc ■  (9,) 

a  ref 

This  is  somewhat  similar  to  the  approach  of  Reference  28.  Equation  (99)  is  used 

with  =  .8  to  establish  at  M^  =  .933.  Equation  (99)  is  then  used  again 

to  establish  C„  at  M  =  1.067  with  M  r  =  1.2.  Interference  corrections  are 
Up  »  ref 

made  for  each  fin  at  M^  =  .8,  .933,  1.06;  (  and  1.2  as  explained  later.  The 
for  the  body  is  added  to  the  lifting  surface  contributions  for  the  above  Mach 
numbers.  The  final  total  value  at  Mw  is  given  by  a  four-point  Lagrange  in¬ 
terpolation. 

2. 2. 2. 1.4  HIGH  SUPERSONIC  FLOW  (M^  >  M£  ) 

*e 

At  higher  Mach  numbers,  where  >  M^  ,  the  potential  model,  for  determining 
the  lifting  surface  roll  damping  is  assumed  to  be  invalid. 


In  analogy  with  the  static  prediction  for  Mw  >  M^  ,  the  strip  loading  is 
assumed  to  be  proportional  to  the  local  chord  and  distance  from  the  centerline, 
lhe  loading  is  forced  to  zero  at  the  tip  in  an  elliptical  manner.  With  these 


assumptions , 
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where 

lx  =  n/4  -  (1  ‘  A)/3, 

12  =  1/3  -  0  -  A)  71/16, 

13  =  71/16  -  2(1  -  \)/15, 

and  C Cjj  )p  is  the  strip  theory  value  for  two  fins.  This  approach  was  evaluated 
in  Reference  4. 

2. 2. 2. 2  PITCH  DAMPING 

The  total  pitch  damping  coefficient,  C  +  Cm(j,  is  made  up  of  two  terms. 
The  first  term,  C  ,  is  due  to  the  rotational  rate  about  the  moment  reference 

mq 

axis.  The  second  term  is  due  to  a  downward  translational  acceleration,  w/V^, 
perpendicular  to  the  body  centerline. 

The  original  method  for  low  Mach  numbers  was  based  upon  small  disturbance 
lifting  theory  (Reference  3).  The  method  developed  by  LMSC  as  presented  in 
References  25  and  26  was  intended  to  apply  over  the  entire  Mach  number  range  at 
a  -  0.  However,  it  is  limited  in  lifting  surface  aspect  ratio  for  low  supersonic 
flows  and  to  one  lifting  surface.  For  this  and  other  reasons,  the  f:hoice  of 
method  for  pitch  damping,  as  earlier  indicated  for  body  pitch  damping,  is  a  pro¬ 
gram  input.  The  LMSC  method  is  based  on  slender  wing-body  methods  in  the  sub¬ 
sonic  and  low  supersonic  flow  regime.  At  high  Mach  numbers,  it  is  based  upon 
Newtonian  embedded  strip  theory. 

Even  if  the  LMSC  option  is  chosen,  the  older  computational  estimate  is 
obtained  first  in  the  program  flow  and  the  newer  estimate  subsequently.  The  LMSC 
model  will  be  presented  in  the  section  on  interference. 

2. 2. 2. 2.1  SUBSONIC  FLOW  (M^  <  .8) 

Rotation  about  the  moment  center  induces  the  following  angle-of-attack  on 
the  wing 

<*(x,y)  =  ^-  [x  +  Cref  (x  -  x  )]  . 

00  ® 


(101) 
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For  the  wing  alone  case,  x  =  0,  =  C,  and  x£g  is  given  C  units.  Otherwise, 

y.  and  x  are  in  caliber  units  and  C  r  =  D. 
eg  ref 

(AC  )  =  (AC  )  q/V  ,  obtained  by  matrix  inversion  for  the  collocation  ap- 

p  q  p  q  00 

proach  of  Equation  (51),  is  an  even-valued  distribution  in  y.  Note  that  the 
second  term  in  Equation  (101)  contributes  an  equivalent  constant  angle-of-attaek 
lift. 

The  pitch  damping  coefficient,  C  ,  is  defined  as 
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Using  the  definition  of  Equation  (102)  one  obtains  for  C 


C  =  -  4 
mq 
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ytan  Yi  +  C(y) 
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s  -  C* 
ref  ref 

The  planform  is  again  extended  through  to  the  body  centerline.  Here,  C ^  is 
assumed  to  be  zero.  A  review  of  Reference  3  indicates  that  various  methods  are 
available  for  estimating  C  £.  However,  the  computational  cost  is  incompatible 
with  the  rapid  estimates  desired  here  and  the  additional  accuracy  unwarranted. 


2. 2. 2. 2. 2  LOW  SUPERSONIC  FLOW  (1.2  <  <  M£  ) 

e 


The  C  problem  proceeds  as  before: 

w'(£,  n)  =  -  ^  U  +  cref  (x  -  xcg)]  . 
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Specific  relations  for  (AC^)^  =  (q/Vw)  (AC^)^  are  given  in  Reference  3  for  the 
variable  part  of  w(£,  r)). 

The  vertical  acceleration  problem  for  a  is  an  unsteady  potential  problem 
which  obeys  the  following  equation  for  the  perturbation  potential: 


»!* .  (H»  - 1)  .  a_  ^ 

3v2  * 


9x2  a2  9t2  a_  3x3t 
00  00 


(105) 


46 


where  is  the  free  stream  value  of  the  speed  of  sound.  The  boundary  condition 


for  the  dimensionless  vertical  velocity  is 
w(£,  H)  =  -  at  . 

Equation  (106)  may  be  rewritten  as 

w(4,  n)  =  -  «[  —  *-  +  t - ] 


(106) 


(107) 


The  first  term  of  Equation  (107)  has  the  form  of  the  variable  part  of  equa¬ 
tion  (106).  The  form  of  Equation  (107)  suggests  a  solution  of  the  form 


/  M2  \  M2 

4>  =  <i  ft  -  —  —  +  a  —  4>2  . 

V  p2  vj  p2 


(108) 


On  the  surface  4>  =  -1,  <t>_  =  -x/V'  .  The  first  term  of  Equation  (108)  satis- 

lz  zz  00 

fies  tne  differential  Equation  (105)  and  leads  to  the  equation 


a2»i  .  p2  ?24»i 


(109) 


Thus  it  is  seen  that  4>!  is  the  lifting  potential  with  a  -  1.  It  is  also  seen 
that  d>2  is  the  potential  corresponding  to  the  variable  part  of  C  with  q  =  1. 
Equation  (108)  may  be  rewritten  as 

/  M2x  \  M2 

4>  =  of  t - W>,  +  a  —  .  (110) 

\  P2vJ  (“=1)  P*  (’=1) 

The  equivalent  of  the  Bernoulli  relationship  for  the  time  dependent  case  is 


V  <b  + 
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(vj  +  fy)  •  (Vj  +  fy) 


+  r-T  =  V2  + 
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Equation  (111),  evaluated  at  t  =  0  from  Equation  (110)  with  the  usual  small 

terms  of  0(^  •  ^<}>)  neglected,  yields  an  estimate  of  (AC  )•  as 

P  W 


M2  4(4)  ,  Mix 

(AC  )•  =  a  -  (AC  )  -  - — - (AC) 
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2. 2. 2. 2. 4  HIGH  SUPERSONIC  FLOW  (M^  >  M£  ) 

e 

A  simple  strip  theory  model  is  considered  here.  Cm(j  is  not  considered. 


8(C  ) 
a  F 

Cmq  =  "  C;efCrb(H-? 


b/2  _  ytan  y  +  C(y) 


«  -  *cg>  Crefl*  dxdy 


0  ytan  y 


^  C2 

2  "(IF  1  ref 


i  /btan  y 
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2. 2. 2. 3  DYNAMIC  DERIVATIVES  FOR  TOTAL  CONFIGURATION 

For  roll  and  pitch  damping,  no  interference  is  assumed  for  the  high  Mach 
number  range.  The  LMSC  pitch  damping  will  be  presented  separately  here. 

2. 2. 2. 3.1  SUBSONIC  AND  LOW  SUPERSONIC  FLOW 

It  has  already  been  noted  that  the  lifting  surface  contribution  was  obtained 
for  a  planform  consisting  of  an  extension  of  the  leading  and  trailing  edges  to 
the  body  centerline. 

An  additional  lifting  surface  interference  is  based  on  slender  body  theory 
(see  Reference  29).  The  lifting  surface  contribution  in  the  presence  of  the  body 


is  given  in  terms  of  a  span  to  diameter  ratio  for  C 

_  r  /d  \ 

tot  w 


£p‘ 


016) 


Here,  however,  (C^p is  based  on  the  modified  planform.  For  a  complete  con¬ 
figuration, 

C£p  =  +  ^ip^WB  *  ^VtB  •  (117:) 

A  similar  fin-body  interference  term  for  C  +  C  •  is  given  as 
'  mq  mof 


(Cmq  *  ChioVb  _  /§_ 

*  'J).  "  8  W 


mq 


ma  D=0 


(118) 


At  Mach  numbers  below  1.2  an  additional  C  ^  term  due  to  downwash  lag  from 
the  canard  onto  the  tail  (see  Reference  23)  is  given  as 


^Sno^T(V) 


2  (x  „  -  x  )  (x  „  -  x  D)  3CVT»w,,N 
v  cpT  cpw/v  cpT  eg  NT(V) 


D2 


(119) 


da 


lor  M„  >  1.2,  (C>i)T(v)  =  0. 

For  a  complete  configuration  CmQ  +  is  given  by 

Cmq  +  Cma  -  ^Cmq  +  Cma^B  +  ^Cmq  +  Cma^WB  *  ^Cmq  + 
+  ^mopT(V)  ’ 


mor  TB 


(120) 


49 


“31 


2. 2. 2. 3. 2  TRANSONIC  FLOW  (.8  <  <  1.2) 

The  approach  for  Cn  and  C  +  C  *  is  the  same.  Values  at  =  .8,  .933, 

SL  p  mq  ma  00 

1.067,  and  1.2  for  the  complete  configuration  are  used  in  a  4-point  Lagrange 

interpolation  to  obtain  total  C.  and  C  +  C  •  at  a  given  M  . 
r  £p  mq  ma  » 


2. 2. 2. 3. 3  LMSC  WING-BODY  PITCH  DAMPING 

For  M  <  1,  slender  wing-body  theory  is  used  where 

(C  +  C  •)_  =  (C  +  C  -)  -  (C  +  C  •)  . 

mq  ma  T  mq  ma  Jl  15  mq  ma  o 

Note  that  this  theory  is  for  a  body-tail  configuration  only 

D,  \2 


where 


(C  + 
v  mq 


C„;>TB  =  -Hr)  0-23  M)j[Sr^-.*3M-.CI] 


-  41  -T 
AD 


(1  -  A  )  -  2* —  tanYi 

rT 


l  nj, 


1  +  2 


n  = 


2At 


3(1  +  \T)  ’ 


x  D  +  C  C  _ 
Xte/D  =  •I_D - -  -  -T  +  C/D  ’ 


K  =2/ 
ma 


AR  =  AR„ 


1  +  AR  p/4  +(P) 


(1  +  DT/bT) 


il  + 


dt  i-atn 


(121) 


(122) 

(123) 

(124) 

(125) 

(126) 

(127) 
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and  AR  is  the  modified  planform  aspect  ratio  obtained  by  extending  the  leading 
and  trailing  edges  tc  the  body  centerline. 


For  supersonic  flow,  1  <  Ma  <^1  +  (ARp/4)2,  Equation  (121)  applies  again. 


/bT\2 

(Cmq  +  CmopTB  =  "  4  \D/  Kma 


te 


D  '  Xcg  -  C/D  l1  *N/Kma 


1  - 


x~  D  +  C  _  \_o 

T  rT  x  ' 2 

- - -  -  Cg 


XT  "  *  « 

T  cg 


(C  +  C  •)„ 
raq  morB 


(128) 


A  maximum  ARp  corresponds  to  ~^1  +  (ARp/4)  2  =  2  or  ARp  <  2.31.  For  larger  aspect 
ratio  planforms  the  code  reverts  to  the  older  theory. 


where 


4*  = 


a* 


flow, 

>1 

00 
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1 
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and  f*  is  the  dynamic  pressure  ratio  across  the  bow  shock  and  is  defined  as 

1.0 


f*  = 


10 


Kn  -  1.25 
12.669724 


.17 


Kn  <  1.25 


1.25  < 


11 


(132) 
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2.3  METHODS  OF  ANALYSIS  SUMMARY 

The  body-alone  static  analysis  methods  are  summarized  in  Figure  13.  The 
Mach  number  region  division  points  are  nominally  at  H#  =  .8,  1.2,  and  Mg  or  M^ 

e 

Methods  for  computing  the  wing-alone  and  interterence  static  aerodynamics 
are  given  in  Figure  14. 
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Finally,  the  methods  to  compute  the  dynamic  derivatives  are  summarized  in 
Figure  15.  Input  option  refers  to  use  of  LMSC  or  other  prediction  methods  for 
C  +  ^ma  ^or  body-alone  or  body-tail  configurations. 
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Figure  15.  Methods  for  Computing  Dynamic  Derivatives 

At  high  angles-of-attack,  the  body-alone  or  body-tail  computations  are 
based  on  the  empiricism  of  Reference  24.  Only  static  derivaties  are  determined. 

Computation  times  depend  on  configuration,  code  option,  and  Mach  number.  A 
computation  time  can  range  from  less  than  a  second  uo  between  30  seconds  and  a 
minute  per  Mach  number  on  the  CDC  6700  computer. 
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3.  APPLICATIONS  TO  CURRENT  DESIGNS 


In  this  section,  data  are  presented  for  a  variety  of  complete  configurations 
of  current  interest.  In  most  cases,  experimental  data  are  shown  for  comparison 
with  the  theoretical  predictions  and  an  estimation  made  of  the  accuracy  of  the 
predictions.  Both  the  static  and  dynamic  coefficients  are  presented  in  many 
cases. 


TACTICAL  WEAPONS  WITHOUT  LIFTING  SURFACES  (BODY  ALONE  AERODYNAMICS] 


Since  the  aerodynamic  coefficients  of  the  individual  body  components  and 
many  of  the  nose^afterbcdy  shapes  are  presented  in  Section  4,  the  data  presented 
here  are  limited  to  complete  configurations.  In  Figures  16  through  20,  the 
static  and  dynamic  aerodynamic  coefficients  are  presented  for  several  spin-sta¬ 
bilized  projectiles.  In  Figute  16,  the  predicted  characteristics  of  the  basic 
Army/Noyy  Spinner  are  shown.  In  general,  the  static  aerodynamics  are  predicted 
reasonably  well  with  the  errors  generally  staying  less  than  10  percent  throughout 
the  Mach  number  range.  The  dynamic  derivative  predictions,  utilizing  the  empiri¬ 
cal  method  (see  Section  2.2.1),  are  generally  poor  and  are  of  the  right  order-of- 
magnitude  only. 

Figures  17  and  18  give  the  estimated  aerodynamic  coefficients  for  a  inch 
Rocket  Assisted  Projectile  (PAP)  and  an  Improved  5  inch/54  Projectile  (Reference 
30),  respectively.  The  5  inch  RAP  has  a  nose  length  of  7.5  calibers,  a  2.2 
caliber  afterbody,  a  .5  caliber  boattail,  and  a  rotating  band  near  the  end  of  the 
afterbody.  The  Improved  round  has  a  2.75  caliber  nose,  a  1.45  caliber  afterbody, 
and  a  1.0  caliber  boattail.  The  Improved  round  employs  a  discarding  rotating 
band  and  therefore  has  no  such  protrusion.’  in  flight.  In  both  of  these  cases, 
there  is  good  agreement  with  experimental  data  obtained  for  the  drag  coefficient 
throughout  the  Mach  number  range  where  data  are  available.  Fair  agreement  was 
achieved  for  the  normal  force  and  center-of-pressure  location.  For  the  5  inch  RAP 
prediction,  the  agreement  approaching  M^  =  2.5  is  very  good,  whereas  at  lower 
Mach  numbers  the  theory  is  consistently  about  10  percent  low.  This  is  due  in 
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Comparison  of  Theory  and  Experiment  for  Army/Navy  Spinner 


Experiment 


part  to  the  inability  of  the  inviscid  theory  to  predict  the  afterbody  lift  cor¬ 
rectly  at  low  supersonic  Mach  numbers.  For  the  Improved  round,  the  dynamic  de¬ 
rivatives'  predictions  are,  in  general,  good.  The  LMSC  methodology  (Reference  25) 
and  the  empirical  predictions  (both  of  which  are  available  in  the  Aeroprediction 
Code  at  the  user's  option)  show  good  agreement  for  C  +  except  in  the  high 
transonic  flow  region,  where  the  agreement  is  fair. 


In  these  cases,  was  generally  set  at  M^  =  2.0  resulting  in  the  use  of 

the  potential  theory  rather  than  the  second-order  shock  expansion  theory  for 

M^  =  2.0.  If  M^  had  been  set  at  M^  =  2.5,  the  potential  theory  would  have  been 

applied  at  this  Mach  number  also,  and  the  predictions  at  M  =  2.5  for  and  CVI 

oo  D  N 

a 

would  have  been  considerably  lower.  This  is  a  result  of  the  Mach  wave  angle 
lying  too  close  to  the  initial  body  slope. 

Figures  19  and  20  present  the  aerodynamic  coefficients  for  two  Army  shells, 
the  175  mm  M437  (Reference  31)  and  the  155  mm  projectiles.  Again,  good  agreement 
is  obtained  for  the  drag  coefficient,  and  fair  agreement  is  achieved  for  the 
normal  force  and  center-of-pressure  location. 


Figures  21  and  22  illustrate  the  Code's  capability  to  predict  the  aerody¬ 
namics  for  blunt  bodies  and  bodies  with  flares.  In  each  case,  it  was  necessary 
to  modify  the  body  geometry  to  be  compatible  with  the  second-order  perturbation 
theory  used  for  low  supersonic  Mach  number  calculations.  The  hemispherical  nose 
on  the  5.5  caliber  hemisphere-cylinder  body  was  modified  to  include  a  0.1  caliber 
conical  frustrum  with  AR  =  .01  at  the  aft  end  (L.,  remained  constant).  The  re- 
mainder  of  the  body  was  described  directly.  The  NASA  flared-body  nose  was  re¬ 
placed  by  a  sohere-cone-f lare  body.  The  spherical  cap  was  allowable  in  this  case 
because  of  the  slight  slope  of  the  afterbody  prior  to  the  flare.  For  the  hemi¬ 
sphere-cylinder,  the  drag  predictions  compare  favorably  with  data  (Reference  32) 
but  the  and  x  predictions  are  "hit  and  miss,"  as  shown  in  Figure  21.  In 
the  case  o?  the  flared  body,  a  slight  increase  in  the  bluntness  resulted  when  the 
nose  was  changed  to  the  spherical  cap.  As  shown  in  Figure  22,  the  effect  of  this 
change  is  particularly  noticeable  in  the  and  predictions.  The  drag  and 
normal  force  are  overpredicted  by  as  much  as  15  percent  in  the  low  supersonic 
region  compared  to  data  (Reference  33).  The  x  predictions  are  within  one-half 
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Figur.  22.  NASA  Flare  Body  Data  Comparison  for  Supersonic  Mach  Numbers 
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small  perturbation  theory  and  the  second-order  shock  expansion  theory.  It  is 
well  known  that  above  M  =  2.0  the  perturbation  theory  begins  to  weaken.  For 
the  blunt  bodies,  the  theory  does  not  hold  at  all.  The  second-order  shock  expan¬ 
sion  theory,  on  the  other  hand,  is  best  applied  in  the  high  supersonic  region. 
In  the  majority  of  cases,  its  accuracy  degrades  rapidly  below  =  3.0.  To  ac¬ 
count  for  this  difficulty  and  to  eliminate  the  discontinuity  in  the  predictions, 
a  weighted  fairing  was  used  in  this  region. 

Projectile  Performance.  In  recent,  years,  there  has  been  a  renewed  emphasis 
on  improving  the  design  of  projectiles.  For  the  A/A  defense  role,  the  major  con¬ 
siderations  are  time  of  flight  and  accuracy.  For  surface-to-surface  applica¬ 
tions,  the  usual  concern  is  maximum  range  as  well  as  accuracy. 

It-  minimizing  the  time  of  flight,  it  is  necessary  to  optimize  the  ballistic 
coefficient,  p  =  w/cj)Srcf  (Reference  34).  The  importance  of  this  parameter  is 
recognized  in  the  following  equations  for  the  velocity  and  time  of  flight  as  a 
function  of  range: 


v  =  vo  «i>  [-it  "*] 

and 


t  = 

f  vopg 


exp 


as 

2B 


(133) 


(134) 


The  drag  coefficient,  C^,  is  assumed  to  be  constant  over  the  range  of  interest. 
For  increasing  values  of  p,  the  velocity  increases  and  the  time  decreases  (both 
exponentially)  for  a  given  range,  R^.  One  would  at  first  feel  that  a  maximum 
value  of  p  would  provide  the  best  performance.  However,  when  designing  a  round 
for  a  given  gun  system,  the  muzzle  velocity  depends  directly  on  the  projectile 
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weight.  Usually  the  new  design  is  replacing  an  old  one  Ii  the  new  weight  is 
less  than  the  previous  standard  weight,  the  muzzle  velocity  for  the  new  design 
can  be  estimated  by  relating  it  to  the  muzzle  energy. 


m  V2  =  m 
o 


std 


V2 

°std 


or 


V 

o 


\  "  /  '  °SLd\  “  J 


(135) 


Equation  (135)  provides  reasonably  good  estimates  as  long  as  the  weight  dif¬ 
ference  is  not  large.  Further  examination  of  Equations  (133)  and  (134)  reveals 
that  increasing  the  projectile  weight  may  increase  P  but  it  also  decreases  V  . 
The  result  is  that  for  a  given  range  and  gun  energy,  an  optimum  value  for  p 
exists  which  will  result  in  the  maximum  velocity,  V,  at  range,  R^,  and  a  minimum 
time  of  flight. 

So  far  the  discussion  has  centered  on  the  projectile  weight.  Obviously,  if 
the  drag  is  minimized,  the  velocity  retardation  is  minimized  and  the  performance 
is  optimized  further.  Therefore,  a  good  approach  to  optimizing  the  projectile 
performance  is  to  first  minimize  the  drag,  then  to  determine  the  optimum  p  for 
the  range  of  interest.  In  the  case  of  intercepting  an  air  target,  the  range 
selected  might  be  the  median  range  over  which  intercept  is  desirable  or  it  may  be 
the  range  at  which  maximum  accuracy  or  advantage  is  required. 

The  above  procedure  was  used  to  optimize  the  performance  of  a  30  mm  projec¬ 
tile.  The  range  of  interest  was  3.0  km.  The  performance  and  physical  character¬ 
istics  of  the  standard  projectile  were  known.  Sample  trajectory  calculations 
indicated  the  median  Mach  number  to  be  about  3.0.  An  L/D  ratio  of  5.0  was  se¬ 
lected  due  to  the  total  round  length  constraints  and  also  to  provide  good  sta¬ 
bility.  Generally,  the  higher  the  fineness  ratio  (L/D),  the  more  difficult  it  is 
to  spin-stabilize  the  body.  With  these  inputs,  the  optimum  shape  for  minimum 
drag  could  be  determined.  The  Aeroprediction  Code  was  then  used  to  obtain  the 
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initial  estimates  of  the  aerodynamic  coefficients  for  the  design.  Figure  23 
shows  the  profile  of  the  body  and  the  predicted  coefficients.  The  optimum  shape 
indicated  has  a  3.2  caliber  nose  and  a  1.8  caliber  boattail  with  a  base  diameter 
of  *7Dre£.  It  should  be  noted  that  the  actual  theoretical  shape  had  a  slight 
curvature  in  the  boattail  near  the  afterbody-boattail  juncture.  Later,  when 
these  projectiles  were  fabricated,  the  profile  was  approximated  by  a  conical 
boattail  for  ease  in  manufacturing.  Comparison  of  the  theory  with  the  ballistic 
range  data  obtained  later  showed  good  agreement  in  each  case.  As  a  result,  the 
predicted  values  could  be  used  with  greater  confidence  throughout  the  Mach  number 
range  to  conduct  a  stability  analysis. 

Having  determined  the  minimum  drag  profile,  it  was  necessary  to  find  the 
optimum  ballistic  coefficient.  Figure  24  shows  the  variation  of  t^.  with  pro¬ 
jectile  weight  for  several  ranges.  The  weight  at  which  t^  is  a  minimum  suggests 
the  optimum  value  for  p.  Note  that  the  optimum  weight  (and  therefore  p)  varies 
with  the  range.  The  variation  in  P  with  the  range  can  be  visualized  in  Fig¬ 
ure  25.  Having  determined  the  desired  weight,  it  is  then  left  to  the  designer  to 
determine  the  fuze  weight  (if  any)  and  the  interior  cavity  which  would  provide 
the  center-of-gravity  location,  polar  (spin  axis)  inertia,  and  the  transverse 
(pitch  axis)  inertia  for  that  weight  which  are  necessary  for  the  stability  analy¬ 
sis.  The  benefits  achieved  through  shape  optimization  can  be  substantial.  In 
our  sample  case,  reductions  in  the  drag  and  time  of  flight  were  34  percent  and  32 
percent,  respectively,  as  compared  to  the  standard  design. 

The  stability  requirements  for  projectile  design  are  discussed  in  detail  in 
in  Reference  34.  Basically,  the  main  concerns  are  the  gyroscopic  stability  factor 
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Figure  24.  Time  of  Flight  of  an  Optimum  Drag  Projectile  vs  Flight 
Weight  to  Slant  Ranges  of  2.0  km,  3.0  km,  and  4.0  km 
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Figure  25.  Minimum  Time-of-Fl ight  Ballistic  Coefficient  for  a 
30  mm  Optimum  Drag  Projectile 
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For  dynamic  flight  stability,  S  >1.0  and  0  <  <  2.0.  In  practice,  good  pro- 

jectile  designs  have,  for  standard  atmosphere  conditions  at  launch,  S  >1.3 

&min 

and  8^  «  1.0.  If  the  values  of  and  are  determined  for  a  particular  design 
to  be  outside  of  these  ranges,  it  is  generally  easier  from  an  aerodynamic  stand¬ 
point  to  look  at  ways  tc  vary  the  I  /I  ratio  and  the  center-of-gravity  location 

x  y 

first  to  correct  the  problem.  If  the  physical  limitations  render  this  approach 
unsuccessful,  then  the  basic  exterior  design  must  be  changed,  usually  through 
increasing  the  length  of  the  cylindrical  section  of  the  body  with  corresponding 
reductions  in  the  boattail  length  or  nose  length.  The  gyroscopic  stability  of  a 
spin-stabilized  projectile  usually  increases  after  muzzle  exit  (at  least  until 
apogee).  Therefore,  it  is  most  critical  that  adequate  stability  be  obtained  near 
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the  muzzle.  However,  both  and  vary  during  the  flight  and  should  be  ex¬ 
amined  throughout  the  trajectory  to  ensure  the  projectile  will  fulfill  its  mis¬ 
sion.  The  values  for  these  parameters  for  the  30  nm  Optimal  Projectile  sample 
case  were  =  2.23  and  =  0.96  for  standard  atmospheric  conditions  at  the 

muzzle.  The  value  for  was  determined  based  on  the  ballistic  range  data  since 

the  theory  predicted  C  +  C  •  and  C  poorly  in  this  case. 

3  r  mq  ma  n  r  3 
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3.2  TACTICAL  WEAPONS  WITH  LIFTING  SURFACES 

Figures  27  through  33  give  the  aerodynamic  coefficients  for  several  missile 
configurations  with  one  or  two  sets  of  lifting  surfaces.  In  each  case,  a  fin 
set  is  comprised  of  four  fins  in  a  cruciform  arrangement. 

In  Figures  26  through  29,  the  aerodynamics  for  configurations  with  only  one 
set  of  lifting  surfaces  are  shown.  The  first,  Figure  26,  is  the  Basic  Finner 
model,  which  has  a  2.84  caliber  conical  nose  and  a  total  length  of  10.0  calibers. 
The  triangular  airfoil,  rectangular  planform  fins  have  AR  =  2.0  and  t/C  =  .08, 
where  C  =  and  AR  =  Aspect  Ratio.  The  geometry  is  representative  of  the  data 

base  and  can  be  modelled  precisely  in  the  Aeroprediction  Code  with  the  exception 
cf  the  blunt  trailing  edge  for  the  fins.  In  the  computer  model,  the  trailing 
edge  was  given  a  radius  of  t/2.  As  expected,  the  theoretical  predictions  match 
the  data  (Reference  35)  quite  well,  including  the  dynamic  derivatives  except  in 
the  transonic  region  where  the  fin  effects  are  underpredicted. 

The  next  three  examples  can  be  classed  as  variations  of  the  Basic  Finner 

♦ 

configuration.  Each  has  a  tangent  ogive  nose,  the  Modified  Basic  Finner  and 
the  Air  Slew  configurations  having  a  2.5  caliber  nose  and  the  TMX-2774  having 
a  3.0  caliber  nose.  The  Modified  Basic  Finner  (Reference  36)  and  Air  Slew 
(AEDC)  have  an  overall  length  of  10.0  calibers  and  a  modified  double-wedge  air¬ 
foil  section  fin  of  AR  =  1.0,  t/C  =  .036,  =  1.33D,  \  =  .5,  and  yj  =  53°.  The 

leading  edge  (LE)  and  trailing  edge  (TE)  bluntnesses  are  small.  The  Air  Slew  is 
a  scaled  version  of  the  Modified  Basic  Finner  except  for  the  fin  bluntnesses. 
For  the  Air  Slew,  the  leading  and  trailing  edges  are  sharp.  The  data  for  the 
Modified  Basic  Finner  at  M  =  2.5  and  the  Air  Slew  at  M_  =  1.3  illustrate  the 
capabilities  of  the  adapted  Martin  High  Alpha  empircial  algorithm.  Noting  the 


small  differences  in  the  configurations,  the  Code  predicts  the  respective  trends 
very  well.  The  slope  and  magnitude  of  the  predictions  compare  quite  well 
except  at  the  higher  Mach  number  and  20°  <  Of  <  40°.  In  this  region,  the  errors 
range  from  5  to  30  percent,  reaching  a  maximum  at  a  =  30°.  The  x  predictions 
show  similar  trends  and  are  of  the  right  order  of  magnitude  only  at  M^  =  2.5. 
The  negative  trend  in  the  prediction  with  increasing  Of  is  not  supported  by 
the  data. 

The  last  of  the  body-tail  configurations  is  the  NASA  TMX-2774  body  (Figure 

29)  which  is  13.06  calibers  long  with  a  .53  caliber,  4.2°  boattail,  and  fins  of 

AR  =  .84,  t/C  -  .07,  Cr  =  1.9,  \  *  .42,  and  Vi  =  62.9°  (Reference  37).  The  LE 

and  TE  are  sharp.  The  data  include  the  effect  of  two  strakes,  180°  apart,  which 

were  not  included  in  the  model  due  to  limits  on  the  number  of  lifting  surfaces 

and  the  small  aspect  ratio.  In  this  case,  the  and  x  predictions  compare 

well.  Again,  it  should  be  pointed  out  that  for  small  aspect  ratios  the  influence 

on  the  opposite  half  wing  tip  is  neglected.  Also,  no  interference  effects  are 

considered  at  high  supersonic  Mach  numbers.  As  a  result,  and  C  +  C  *  pre- 

lip  mq  moi  * 

dictions  for  small  AR  often  behave  erratically.  The  peak  in  the  curve  for  C„  is 
representative  of  the  problem  and  can  be  seen  elsewhere. 

Figures  30  through  32  present  dat3  for  configurations  having  two  sets  of 
lifting  surfaces.  The  large-winged  (AR  =  1.52)  TMX-1751  conf iguration  (Fig¬ 
ure  30)  incorporated  body  strakes  which  were  included  in  the  data  but,  as  before, 
were  not  included  in  the  model  (Reference  37).  This  partially  accounts  for  the 
differences  in  and  x  .  The  agreement  with  data  for  the  TMX-3070  (Figure  31) 

is  very  good  foraCp  and  in  spite  of  protrusions  on  the  actual  model,  includ¬ 
ing  a  fin  assembly  sleevea(Reference  37)  The  prediction  for  x  is  fair  to 
poor,  however,  with  errors  ranging  from  .15  calibers  (1.0  percent  of  the  distance 
from  the  nose)  to  4.9  calibers  (35  percent).  This  is  probably  caused  by  an  un- 
detprediction  of  the  lift  for  the  highly  swept  canards  and  errors  in  the  pre¬ 
dicted  pressure  distributions  on  the  hemispherical  nose  noted  earlier.  Note 
again  the  fluctuation  in  due  tc  the  lower  .AR  tail.  The  next  example,  in 
Figure  32,  is  an  illustration  of  probably  the  most  complex  configuration  which 
can  be  handled  by  the  Aeroprediction  Code.  The  body  i  comprised  cf  a  uph^rical- 
ly  blunted  2.1  caltbjr  nose,  a  cylindrical  afterbody,  and  an  8.9°  boattail.  The 
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nose  is  divided  into  two  distinct  ogival  nose  sections.  There  are  two  sets  of 
swept  lifting  surfaces,  each  having  different  fin  geometries  and  deflections. 
The  Aeroprediction  Code  provides  data  for  the  individual  body  components.  The 
data  shown  in  the  example  are  for  the  total  configuration  with  a  5°  canard  de¬ 
flection. 

One  capability  of  the  Code  which  has  not  yet  been  mentioned  is  illustrated 
in  the  following  example.  In  this  case,  the  Aeroprediction  Code  was  modified  to 
compute  the  total  sectional  normal  force  coefficients  along  a  missile  body  as  a 
function  of  Mach  number  and  angle-of-attack.  The  results  of  the  aerodynamic 
analysis  were  then  used  in  a  finite  element  structural  analysis  to  predict  struc¬ 
tural  integrity.  The  results  of  the  analyses,  shown  in  Figure  33,  indicate  that 
structural  failure  of  the  two  configurations  examined  can  occur  during  maneuver¬ 
ing  flight  due  to  the  large  aerodynamic  bending  moment.  Unpublished  flight  test 
results  confirmed  the  predictions. 
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Figure  33.  Structural  Integrity  Study  Using 
the  Aeroprediction  Code 
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Missile  rmance.  Designing  a  missile  for  optimum  performance  is  a  very 
complex  procedur;  with  a  number  of  tradeoffs  and  considerations.  A  few  factors 
which  are  considered  include  the  velocity,  range,  time  of  flight,  altitude, 
launch  angle,  payload,  maneuverability,  responsiveness,  lift-to-drag  ratio,  and 
performance  of  the  propulsion  system.  In  contrast  to  the  projectile  performance 
optimization  where  the  ballistic  coefficient  plays  an  important  role,  a  missile 
design  is  typically  optimized  for  one  or  more  of  the  factors  listed.  Moreover, 
there  are  more  options  available  to  the  aerodynamicist  in  achieving  the  perform¬ 
ance  goals  due  to  the  lifting  surfaces.  The  Aeroprediction  Code  is  ideal  to 
parametrically  examine  the  aerodynamic  factors  affecting  performance  as  a  result 
of  its  relative  low  cost  and  accuracy. 

To  illustrate,  a  typical  design  approach  is  outlined  which  can  be  used  to 
arrive  at  a  preliminary  surface  launched  missile  design  which  satisfies  the  de¬ 
sired  performance  objectives. 

Initially  the  missile  performance  requirements  must  be  established.  The 
specified  requirements  should  include  maximum  and  minimum  range,  missile 
velocity,  altitude  ceiling,  maneuverability  (g's)  required,  and  airframe  time 
constant.  These  requirements  will  usually  be  driven  by  the  threat  spectrum  to 
be  engaged.  The  launch  and  payload  weight  and  volume  constraints  may  also  be 
specified  due  to  the  given  launching  tube  or  canister  envelope.  The  maximum 
range  requirement  can  be  driven  by  several  factors  including  the  fire  control 
radar  detection  range  and  the  tiring  doctrine  to  be  utilized.  The  minimum  range 
(must-kill  range)  is  usually  established  by  doctrine.  In  this  example,  the 
missile  is  assumed  to  be  launched  from  a  square  canister  of  specified  length,  L. 
The  payload  is  assumed  to  include  the  guidance  electronics,  radome,  warhead,  and 
fuze.  The  average  missile  velocity  is  assumed  to  be  in  the  Mw  =  3.0  range. 

The  maximum  flight  range  for  the  missile  is  a  function  of  the  time  history 
of  the  missile  velocity  which  in  turn  is  the  integral  of  the  thrust  and  drag 
along  the  flight  path. 


The  configuration  sizing,  therefore,  involves  design  tradeoffs  of  the  air¬ 
frame  drag  and  weight  and  the  propulsion  system  thrust  required  to  deliver  the 
payload  to  the  maximum  range  with  adequate  intercept  velocity  and  maneuver  capa¬ 
bility.  Since  additional  missile  weight  adversely  impacts  the  propulsion  system, 
flight  performance,  and  packaging  requirements,  the  designer  generally  desires  to 
determine  the  minimum  launch  weight  to  meet  stated  requirements. 


The  missile  intercept  speed  affects  several  aspects  of  guidance,  including: 
the  maximum  allowable  target  crossing  angle,  which  is  a  function  of  the  target 
speed;  the  seekers  maximum  look  angle  and  the  airframe  angle-of-attack;  the  miss 
distance,  which  tends  to  increase  as  the  closing  velocity  increases  due  to  re¬ 
ceiver  and  seeker  servo  noise;  and  the  missile  guidance  time  constant,  T  .  The 
total  missile  guidance  time  constant  is  defined  as  the  time  from  the  initial 
tracking  loop  input  upon  target  acquisition  to  the  time  it  takes  the  missile  to 
achieve  approximately  63  percent  of  the  desired  rate-of- change  of  the  flight  path 
angle,  Y-  The  missile  guidance  time  constant  can  be  approximated  to  first  order 
using  the  relation 


t 

g 


\  +  'f  +  taf  +  N 


(138) 


where 


Tj  is  the  tracking  loop  time  constant, 

is  the  guidance  filter  time  log, 
r 

is  the  airframe  time  constant, 

N  is  the  effective  navigation  ratio, 

Vc  is  the  missile  closing  velocity, 

Vm  is  the  missile  average  velocity, 

a/Y  is  the  turning  rate  time  constant,  and 

R^e  is  the  radome  boresight  error  slope. 
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A  low  value  of  T  is  desirable  in  order  to  reduce  miss  distance  sensitiv- 
g 

ity  to  target  maneuver  and  provide  adequate  terminal  accuracy  under  adverse  con¬ 
ditions.  The  adverse  conditions  such  as  electronic  countermeasures  or  weather 
result  in  a  reduction  in  the  range  available  for  guidance.  A  useful  criterion 

for  seeker  end  game  performance  is  to  allow  a  minimum  of  10  I  (10  guidance  time 

g 

constants)  to  correct  terminal  heading  error.  It  can  be  seen  from  the  above  re- 

Vc  a 

lation  that  decreases  in  —  and  both  of  which  decrease  with  increasing  missile 

vra  i 

velocity,  decrease  the  total  time  constant.  Also,  limitations  due  to  the  battle- 
space,  which  is  the  distance  between  the  maximum  and  minimum  ranges,  and  the  fir¬ 
ing  doctrine  (such  as  a  shoot-look-shoot  doctrine)  encourage  the  selection  of  as 
high  a  missile  velocity  as  possible.  However,  the  increase  in  the  miss  distance 
as  the  closing  velocity  increases  forces  a  compromise  on  the  selection  of  the 
desired  missile  velocity.  In  our  example,  an  Mw  =  3  missile  velocity  appears  to 
meet  these  restraints. 


The  sizing  of  the  missile  is  a  function  of  the  propulsion  type  and  capa¬ 
bility,  the  type  of  aerodynamic  control  (or  thrust  control  if  used),  and  the 
maneuverability  (g's)  required.  The  propulsion  system  selection  is  based  on 
altitude,  thrust,  and  burn  time  requirements.  Both  simple  boost  and  boost- 
sustain  motors  should  be  examined  and  correlated  with  the  missile  weight  to  pro¬ 
vide  the  required  average  or  sustain  velocity  over  the  desired  range. 

The  choice  of  the  type  of  aerodynamic  control  is  largely  a  tradeoff  involv¬ 
ing  the  location  of  the  controlling  surfaces  on  the  body.  Canards  offer  advan¬ 
tages  in  obtaining  a  quick  response  from  the  missile.  However,  canards  are 
angle-of-attack  limited  and  often  lead  to  pitch  and  roll  control  problems  due  to 
the  interference  of  the  shed  vortices  on  the  rear  stabilizing  surfaces.  Wings, 
which  must  be  Larger  to  produce  the  same  moment  due  to  their  closeness  to  the 
center  of  gravity,  have  similar  advantages  and  disadvantages  as  canards.  Tail 
alone  control  has  the  advantages  of  lower  overall  weight  and  better  piteh/roll 
control  due  to  the  elimination  of  vortex  interaction  with  vortices  shed  from 
lifting  surfaces  located  forward  of  the  tail.  It  has  the  disadvantages  of  longer 
missile  response  time  because  the  control  force  to  rotate  the  missile  is  opposite 
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to  the  desired  direction  of  travel  and  actuator  packaging  must  be  accomplished 
around  the  rocket  motor.  In  this  case,  tail  control  is  selected  for  the  concep¬ 
tual  missile  design.  A  jettisonable  set  of  jet  vane  controls  may  be  added  in 
order  to  implement  the  quick  turn  maneuver  from  the  vertical  when  a  vertical 
launch  system  is  utilized. 

Having  made  the  above  choices,  the  missile  diameter  and  fineness  ratio 
can  be  determined.  Using  an  average  overall  packing  density  of  .062  lb/in3, 
relationships  between  the  missile  length,  diameter,  and  weight  for  the  body- 
tail  type  configuration  can  be  obtained  as  shown  in  Figure  34.  Based  on  these 
curves,  an  8-inch  diameter  configuration  apparently  would  best  meet  the  length 
and  weight  constraints  shown. 


WL,  LAUNCH  WEIGHT  -  lbs 

Figure  34.  Missile  Physical  Characteristics 

Sufficient  missile  maneuverability  is  required  in  the  terminal  phase  of 
flight  to  compensate  for  guidance  errors,  heading  error,  and  target  maneuvers. 
An  estimate  of  the  missile  maneuverability,  am,  required  to  correct  for  guidance 
errors  due  to  noise  and  target  maneuvers  is  provided  by  the  relation 

am  >  (3aT  +  8)  g  (139) 
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where  a^,  is  the  maneuvering  target  acceleration. 

The  miss  distance  due  to  initial  heading  error,  KE0,  is  based  on  an  average 
missile  acceleration,  a^,  applied  over  time,  t^,  and  is  given  by 

"HE  =  Vo  <>  '  V t  <14« 

o 

where  t  »  10  t  and  the  initial  missile  acceleration  (assuming  no  g  saturation) 
o  g 

is : 


a 

m 


N  V 

Ttf  ■*.  <8's>' 
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For  a  2.0  .second  guidance  time,  T  =  .2,  and  N  =  3,  an  initial  HE  of  10°  can 

8 

be  corrected  and  the  miss  distance  eliminated  provide!  the  missile  can  achieve  a 
27  g  maneuver. 

Figure  35  shows  the  radome  boresight  error  slope  variation  versus  nose  drag 
coefficient  for  several  nose  shapes  and  fineness  ratios.  To  minimize  the  bore- 
sight  error  while  achieving  relatively  low  drag,  a  2.25  caliber  Von  Karman  nose 
contour  is  selected. 


It  is  now  possible  to  define  the  overall  missile  fineness  ratio.  For  the 
trimmed  condition,  it  is  necessary  that  the  summation  of  the  moments  be  zero 
such  that 


Cmo<  aTrim  +  Cin^  ^F 


(142a) 


or 


(142b) 


Therefore,  to  obtain  maximum  maneuverability,  it  is  desired  to  minimize  the 
total  untrimmed  pitching  moment  coefficient,  CmQ.  Similarly,'  the  normal  force 
for  the  trimmed  condition  is 


.  =  CVI  a_  .  +  CVI 

Trim  N- 
Trim  a  o 


(142c) 
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Figure  35.  Forebody  Drag  Coefficient  vs  Boresight  Error  Slope 
for  Various  Shapes  and  Fineness  Ratios 

Since  deflection  of  the  tail  surfaces  produces  a  negative  normal  force  tj 
achieve  the  desired  angle-of-attack,  it  can  be  seen  that  the  body  untrimmed 
required  is  greater  than  the  trimmed  condition.  With  this  in  mind,  several  tail 
sizes  and  planforms  were  examined  to  obtain  a  first  approximation  of  the  tail 
size  and  to  estimate  the  negative  normal  force  effect.  Figure  36  shows  the  vari¬ 
ation  of  the  tail  normal  force  due  to  deflection,  C^,  and  the  total  configura¬ 
tion  pitching  moment,  Cma,  for  a  hypothetical  body  versus  taper  ratio  assuming  a 
diamond  airfoil  cross  section.  The  optimum  tail  for  maneuverabililty  is  that 
which  produces  the  least  negative  force  with  deflection  and  the  minimum  total 

pitching  moment.  For  this  example,  the  optimum  tail  was  that  having  a  C./C  = 

-1  t  r 

.364  resulting  in  a  =  -0.67  (deg  ,  based  on  the  planform  area,  S^,) .  If 

the  center  of  gravity  of  the  total  missile  configuration  after  burnout  is  assumed 
to  be  .54L  from  the  nose  tip,  the  normal  force  required  versus  angle-of-attack 
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for  configurations  of  varying  L/D  ratios  can  be  determined.  In  this  case,  the 
total  trimmed  required  for  the  required  maneuver  is  3.4  for  a  maximum  angle- 
of-attack  of  18°.  This  indicates  that  a  body  alone  normal  force  coefficient  of 
4.0  must  be  obtained.  This  is  possible  with  a  L/D  =  18  which  was  shown  earlier 
to  also  meet  the  length  and  weight  constraints. 

The  resulting  preliminary  missile  configuration  and  the  estimated  aerody¬ 
namic  coefficients  are  shown  in  Figure  37  and  38.  Using  these  data  and  the  pro¬ 
pulsion  data,  the  missile  performance  can  be  predicted  for  various  conditions. 
A  typical  trajectory  is  shown  in  Figure  39. 


Figure  39.  Example  Missile  Configuration  Performance 


4.  AIDS  TO  DESIGN 


As  can  be  readily  seen  in  the  material  of  Section  2,  a  large  number  of 
methods  and  data  have  been  brought  together  into  one  computer  program.  The 
resulting  configurational  generality  possible  allows  one  to  examine  the  effects 
of  design  modifications  of  various  configuration  components  on  the  total  con¬ 
figuration  aerodynamics  and  stability.  Handbooks  have  been  previously  published 
(References  5,  22,  and  32)  which  have  provided  aids  to  the  missile  designer  for 
certain  classes  of  configurations.  In  this  section,  a  number  of  design  charts 
are  presented  which  were  generated  using  the  Aeroprediction  Code  for  various  con¬ 
figuration  components  of  both  pointed  and  blunt  bodies.  These  charts  represent 
portions,  reproductions,  and  extensions  of  the  charts  presented  in  these  earlier 
handbooks.  The  user  is  encouraged  to  compare  the  data  presented  with  that  of 
the  earlier  works  where  possible.  The  basic  desire  of  the  authors  is  to  show 
the  trends  due  to  modifications  in  the  design  components  which  will  guide  the 
designer  in  the  preliminary  design  process.  For  body  component  parameters  within 
the  limits  of  the  charts  presented  here,  the  reader  can  use  the  method  of  com¬ 
ponent  superposition  (i.e.,  to  add  up  the  various  effects  of  each  component)  to 
obtain  the  aerodynamics  of  the  complete  configuration. 

4.1  BODY  ALONE  DESIGN 

The  aerodynamic  coefficients  were  determined  for  a  large  number  of  body 
alone,  axisymmetric  shapes  throughout  the  Mach  number  range.  The  characteristics 
of  the  body  components  were  varied  systematically  to  establish  the  trends  in  the 
aerodynamic  coefficients.  Particular  attention  was  given  to  bodies  with  conical, 
tangent-ogive,  and  secant-ogive  nose  contours  as  shown  in  Figure  40  because  of 
their  primary  interest  to  missile  designers. 

4.1.1  DRAG 

As  indicated  in  Section  2,  the  total  drag  of  the  body  is  the  sum  of  the  nose 
wave  drag,  the  skin-friction  drag,  the  boattail  wave  drag,  the  base  drag,  and  the 
drag  of  protrusions,  plus  the  crossflow  contribution.  In  this  section,  data  are 
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Figure  40.  Typical  Nose  Contours  (LN  =3.0  cal) 

presented  for  a  0°.  For  moderate  angles-of-attack,  a  crossf low  drag  correction 
is  presented. 

The  nose  wave  drag  is  influenced  primarily  by  nose  shape,  fineness  ratio 
(L^/D) ,  and  Mach  number.  The  effect  of  these  parameters  on  the  wave  drag  in 
supersonic  flow  is  illustrated  in  Figure  41  for  pointed  bodies.  In  general, 
the  wave  drag  decreases  for  increasing  nose  fineness  ratio  and  Mach  number.  Of 
the  nose  contours  shown,  the  secant  ogive  is  shown  to  provide  the  lowest  wave 
drag  within  the  Mach  number  range  1.0  <  M  <  6.0.  Earlier  studies  (References  5 
and  38)  have  shown  that  the  minimum  wave  drag  for  secant  ogives  for  low  super¬ 
sonic  Mach  numbers  (M  <  2.5)  is  obtained  by  using  a  radius  of  curvature  twice 
that  used  for  a  tangent  ogive  of  the  same  length.  In  this  comparison,  all  of  the 
secant  ogive  profiles  had  a  constant  shoulder  angle  of  6.0°.  For  low  supersonic 
Mach  numbers  (M^  <  2.0)  slender  body  theory  predicts  the  minimum  drag  profile  to 
be  the  well-known  Von  Karman  ogive.  Newtonian  theory  predicts  that  the  conical 
nose  shape  would  become  the  minimum  drag  profile  as  M^  -*  «.  Therefore,  in 
general,  the  expected  trend  in  supersonic  flow  is  that  for  a  given  nose  length 
and  diameter  the  minimum  drag  profile  would  vary  slowly  from  the  Von  Karman 
ogive  to  the  secant  ogive  (with  radius  twice  that  of  a  tangent  ogive)  and  then 
to  a  conical  nose  shape  as  the  Mach  number  increases.  Although  not  shown,  the 
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Figure  41.  Wave  Drag  Comparison  for  Pointed  Bodies 


power-law  nose  contours  also  closely  follow  these  trends  and  may  in  some  in¬ 
stances  exhibit  lower  drag  than  the  secant  ogive  due  to  nose  bluntness.  For 
example,  the  6.0°  (shoulder  angle)  secant  ogive  can  be  approximated  by  the  two- 
thirds  power-law  nose.  As  the  Mach  number  increases,  the  power  used  to  obtain 
the  low  drag  shape  increases  to  a  value  of  1  which  corresponds  to  the  conical 
shape . 

For  a  nose  of  fineness  ratio  of  3.0,  the  variation  in  the  wave  drag  with 
nose  shape  is  shown  in  Figure  42.  In  this  case,  the  more  blunted  tangent  ogive 
is  shown  to  have  the  lower  wave  drag  for  M^  =  3.0.  For  subonic  Mach  numbers,  the 
wave  drag  contribution  is  assumed  to  be  zero  for  M^  <  .5.  The  exact  point  at 
which  the  wave  drag  disappears  depends  on  the  nose  fineness  ratio.  This  is  il¬ 
lustrated  in  Figure  43  which  shows  the  transonic  wave  drag  for  various  tangent 
ogives . 

The  nose  wave  drag  is  also  influenced  by  bluntness  as  shown  in  Figures  44 
and  45  for  cones.  The  data  presented  are  for  a  7.125°  serai-angle  truncated  cone 
with  fineness  ratio  4.0  when  not  blunted  and  for  a  10°  semi-angle  spherically 
capped  cone.  When  adding  the  spherical  cap  or  truncating  the  nose,  the  nose 
length  was  reduced  back  to  the  point  where  the  proper  value  for  (the  radius 
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Figure  42.  Wave  Drag  Comparison  for  Pointed 
Bodies  (Nose  Length  =3.0  cal) 


Figure  43.  Transonic  Wave  Drag  of  Tangent  Ogives 


.6 


of  the  cap  for  the  spherically  capped  cones)  would  be  obtained  while  keeping 
the  slope  of  the  conical  frustrum  constant.  As  can  be  seen,  truncation  alone 
has  a  slightly  greater  effect  than  truncation  with  a  spherical  cap.  Note, 
however,  that  a  portion  of  the  difference  can  be  attributed  to  the  slightly 
different  nose  lengths  obtained  in  the  two  cases  for  the  same  bluntness  ratio. 
The  diffr  >  it  nose  length  results  from  both  the  different  slopes  and  the  smaller 
be  '/  r  'Hus  •  t  the  point  of  tangency  of  the  spherical  cap.  In  either  case, 
bl  11  ....  the  nose  up  to  Rjj/^re£  =  *1  has  a  negligible  effect  on  the  wave  drag 
(or  on  the  other  aerodynamic  coefficients,  as  shown  in  the  following  subsec¬ 
tions).  Further  increases  in  the  bluntness  ratio  results  in  pronounced  increases 
in  th.1  •’<  .g 

Ja  v.-ry  low  fineness  ratio  pointed  noses  with  little  or  no  curvature,  the 
high  body  slope  row  result  in  a  viscous  separation  drag  at  subsonic  and  transonic 
Mach  numbe  n.  This  drag  results  from  a  separation  bubble  formed  by  flow  detach¬ 
ment  at  th.  nose  tip  due  to  strong  adverse  pressure  gradients  and  reattachment 
downstream.  The  phenomenon  is  somewhat  Mach  number  dependent  and  is  clearly 
noticeable  Ccr  Ma  <  .8.  Its  presence  is  apparent  up  to  M^  ~  .95  or  so  depending 
on  the  lone  ?t.gle.  Figure  46  shows  the  magnitude  of  the  nose  separation  drag 
for  coues  with  and  without  afterbodies  at  subsonic  speeds.  The  effect  of  the 
nose  separation  bubble  is  noticeably  reduced  by  the  influence  of  the  afterbody. 
The  effect  of  nose  curvature  is  to  eliminate  the  adverse  pressure  gradients 
and  therefore  the  phenomenon  does  net  occur. 

For  angles-of-attack  near  zero,  the  contribution  of  cylindrical  afterbodies 
to  the  wave  drag  is  assumed  to  be  zero  and  thus  no  additional  contribution  to 
the  wave  drag  occurs.  Flares  and  boattails  contribute  positively  to  the  wave 
drag.  Wave  drag  data  for  flares  are  presented  in  References  5  and  22,  and  the 
reader  is  referred  to  these  publications  in  these  instances. 

The  wave  drag  of  the  boattail  is  primarily  influenced  by  the  boattail  slope 
or  angle,  the  boattail  length,  Mach  number,  and  the  length  of  the  preceding 
afterbody  as  this  affects  the  degree  of  pressure  recovery  on  the  afterbody  and 
subsequent  expansion  at  the  afterbody-boattail  juncture.  For  short  afterbodies 
(L  <  4.0  calibers),  the  pressure  recovery  and  thus  the  boattail  wave  drag  are 
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Figure  46.  Separation  Pressure  Drag  on  a  Cone-Cylinder 

also  influenced  by  the  nose  shape  (Reference  14),  the  influence  generally  becom¬ 
ing  stronger  as  the  nose  shoulder  angle  increases.  In  this  study,  the  effects  of 
nose  shape  and  afterbody  length  on  the  boattail  were  not  examined.  However,  it 
should  be  noted  that  the  numerical  marching  technique  used  in  the  Aeroprediction 
Code  would  allow  one  to  examine  these  parameters,  at  least  to  the  end  of  the 
afterbody. 

The  wave  drag  coefficient  for  conical  boattails  in  supersonic  flow  is  given 
in  Figure  47.  For  a  given  boattail  length,  decreasing  the  base  diameter  in¬ 
creases  the  boattail  wave  drag.  However,  as  can  be  seen  from  Equation  (136), 
decreasing  the  base  diameter  also  results  in  a  lower  base  pressure.  Therefore, 
an  optimum  configuration  with  respect  to  drag  would  balance  these  drag  contri¬ 
butions.  The  overall  decrease  in  total  drag  of  a  configuration  can  be  strongly 
influenced  by  the  boattail  design.  This  is  illustrated  clearly  in  Figure  48 
which  shows  the  variation  in  the  total  drag  of  a  10.5  caliber  body  having  a  3.0 
caliber  tangent-ogive  nose,  a  6.0  caliber  cylindrical  afterbody,  and  a  1.5  cali¬ 
ber  conical  boattail  with  various  base  diameters.  In  subsonic  flow,  decreasing 
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Figure  47.  Wave  Drag  Coefficient  of  Conical 
Boattails  at  Supersonic  Speeds 


base  diameter  continues  to  decrease  the  total  drag.  In  transonic  and  supersonic 
flows,  however,  decreasing  the  base  diameter  to  a  ratio  of  =  .4  or  lower 
results  in  an  increase  in  the  total  drag.  In  supersonic  flow,  as  the  Mach  number 
increases,  the  influence  of  the  base  drag  decreases  rapidly,  thus  allowing  larger 
base  diameters  to  achieve  a  minimum  drag  boattail  design. 

Many  projectiles  also  retain  rotating  bands  in  flight.  In  the  absence  of 
fraying,  the  drag  increment  due  to  a  rotating  band  has  been  functionalized  (Ref¬ 
erence  1)  and  can  be  obtained  by  using  Equation  (17)  and  the  curve  shown  in  Fig¬ 
ure  49.  More  recently,  rotating  bands  have  been  designed  for  full-bore  pro¬ 
jectiles  which  discard  upon  muzzle  exit,  thus  eliminating  this  drag  contribution 
and  allowing  for  a  cleaner  shape. 
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Figure  49.  Rotating  Band  Drag  -  C. 

hRB 

For  a  >  0,  additional  drag  results  from  the  viscous  crossflow  velocity 
(V  =  V  sin  a).  For  moderate  angles-of-attack,  one  needs  first  to  determine 
the  total  normal  force  (as  outlined  in  the  following  section  beginning  on  page 
98).  Then,  the  total  drag  is 

(L  =  C.  cos  ot  +  CL  sin  a  (143) 

DA  N 

where 

C.  =  CL  cos2a.  (144) 

A  D 

o 

The  skin  friction  component  of  the  total  drag  is  computed  using  Van  Driest 's 
method  as  outlined  in  Section  2.  The  skin  fricton  drag  is  influenced  primarily 
by  the  Reynolds  number,  the  heat  transfer,  the  Mach  number,  and  the  point  on  the 
body  at  which  the  boundary  layer  transitions  from  laminar  to  turbulent  flow.  For 
the  body  alone,  the  boundary  layer  will  generally  be  turbulent  over  roughly 
90  percent  of  the  body,  depending  on  flow  conditions  and  overall  body  length. 
For  very  short  bodies  and  for  fins,  the  percentage  will  likely  be  much  less. 
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The  calculation  of  the  skin  friction  coefficient  using  Van  Driest's  method 
is  somewhat  cumbersome.  A  rough  approximation  can  be  obtained  by  using  the 
simple  formulation  (from  Reference  5): 


w 


F  S 


ref 


(145) 


where 


S  is  the  wetted  surface  area, 
w 

S  x  is  the  reference  arc,  and 
ref  * 

C„  is  the  skin  friction  coefficient  determined  from  Figure  50  multiplied  by 

I? 

the  factor  1.15  for  a  body  of  revolution. 


No  attempt  is  made  here  to  examine  the  trends  for  the  base  drag  or  the 
various  parameters  affecting  base  drag  in  detail.  Deviations  in  the  base  drag 
could  occur  with  variations  in  the  nose  shape,  afterbody  length,  angle-of-attack, 
Reynolds  number  (if  sufficiently  low),  surface  temperature,  and  boattail  shape. 
However,  the  most  common  method  to  determine  the  base  drag  in  the  absence  of  any 
mass  base  flow  (such  as  base  bleed  or  rocket  exhaust)  is  an  empirical  method 
based  on  data  for  a  long  cylindrical  afterbody  and  a  fully  turbulent  boundary 
layer  ahead  of  the  base  (Reference  1)  and  adjusting  for  the  presence  of  a  boat- 
tail.  The  relation  for  determining  the  base  drag  in  this  situation  [using  Equa¬ 
tions  (13)  and  (20)  from  Section  2)]  is 


<VRref>3 


(18) 


where  Cp  is  determined  from  Figure  51. 

BA 

For  small  angles-of-attack  and  low  supersonic  Mach  numbers,  the  relation 

AO aba  =  (-.64930  -  .002833MW) (Rg/Rref)3  (20) 

can  be  used  to  adjust  the  base  axial  force. 
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FRICTION  COEFFICIEN 


4.1.2  NORMAL  FORCE  AND  CENTER  OF  PRESSURE 


At  zero  angle-of-attack,  an  axisymmetric  body  exhibits  no  normal  force. 
However,  with  even  small  increases  in  the  angle-of-attack,  the  normal  force 
increases  rapidly.  The  slope  of  the  normal  force  curve  (C^  versus  a),  expressed 
as  the  coefficient  and  its  location  of  action  (center  of  pressure,  x  )  rela¬ 
tive  to  the  body  center  of  gravity  is  of  primary  concern  in  assessing  the  stabil¬ 
ity  of  a  flight  vehicle.  Projectiles  without  lifting  surfaces  are  generally 
statically  unstable  and  must  be  stabilized  through  gyroscopic  means.  For  a 
statically  stable  body,  the  x  must  be  rearward  of  the  x£  .  This  is  often 
achieved  by  placing  lifting  surfaces  on  the  rearward  portion  of  the  body.  The 


location  is  more  or  less  positioned  depending  on  flight  conditions  and  the 


degree  of  manueverability  desired. 


The  normal  force  derivative,  C„  .  and  the  x  for  various  nose-afterbody 

*  cp 

combinations  at  zero  angle-of-attack  are  given  in  Figures  52  through  57.  The 

data  presented  are  for  bodies  with  tangent-ogive,  secant-ogive,  and  conical 

nose  shapes.  The  effect  of  nose  length,  afterbody  length,  and  Mach  number  can  be 

readily  assessed  from  these  graphs.  For  p/L^  >  .9  (high  oupersonic  flow),  the 

trends  for  the  values  of  and  x  are  quite  consistent.  Generally,  for  a 

given  I^j  and  Mach  number,  increasing  afterbody  length  results  in  a  corresponding 

increase  in  C„  and  a  decrease  in  the  x  distance  from  the  nose  (in  terms  of  the 

\  cp 

percentage  of  the  total  body  length).  Note  that  the  data  presented  are  for 
a  =  0°.  As  a  result,  for  a  tangent  ogive-cylinder  body  and,  to  a  lesser  degree, 
for  the  secant  ogive-cylinder  body,  the  data  indicate  that  on  increase  in  the 


afterbody  length  beyond  L^/ =  3.0  does  not  significantly  increase  .  This 
would  not  necessarily  be  the  case  for  higher  angles-of-attack.  The  effects  of 
afterbody  length  and  Mach  number  are  more  easily  visualized  in  Figures  58 
and  59  for  =  4.0  caliber  tangent  ogive-cylinder  and  cone-cylinder  bodies. 
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Figure  52.  Supersonic  Normal  Force  Coefficient  Variation 

for  Cone  Cylinders 
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Figure  55.  Supersonic  X  Variation  for  Secant-Ogive  Cylinders 
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For  increasing  nose  length  or  decreasing  Mach  number  (0/L^  decreasing), 

increases  gradually  to  a  maximum  before  falling  off. 
a 

In  the  region  0/Ljj  >  .9,  the  results  shown  compare  favorably  with  those  of 
Reference  22.  In  the  transonic  and  subsonic  regions,  the  0/L^  correlation  does 
not  continue  to  hold.  (This  portion  of  the  graph  is  shown  for  illustrative 
purposes  only  and  should  not  be  used  for  data  extraction.)  Note  that  the  pre¬ 
dicted  trend  for  ,  although  not  well  defined,  is  for  Cjj  to  decrease  in  each 
case,  as  0/L^  decreases  below  a  value  of  1.0,  to  a  value  of  2.2,  which  is  about 
10.0  percent  above  the  slender  body  value. 


For  subsonic  flow,  the  total  inviscid  normal  force  for  a  nose-cylinder  body 
is  simply 


CN  =  <CS  >X  +  <CN  >A  (146) 

a  a  a 

where  the  nose  contribution  is 

(CN  =  Cl  tan  6 *  +  Cz  047) 

and  (C^  )^  is  the  afterbody  contribution  (Reference  1).  '^e  values  for  Cj  and 

C2  are  given  in  Figure  60  and  (C^  )A  is  found  from  Figure  61.  The  angle  6*  is 

the  nose  shoulder  angle  at  the  nose-afterbody  juncture.  Slender  body  theory 

yields  the  x  location  for  the  nose  in  subsonic  flow  as 
cp 


(x  ).,  =  Lv,  - 
cp  N  N 


(v°i)» 


ref 


(148) 


For  the  afterbody,  it  is  suggested  that  one  use  Figure  62  taken  from  Reference  1. 
A  simplified  expression  for  the  x  for  the  total  body  is  then 


Xcp  =  ^cp^N  (CN  +  *xcp^NA  (CN  V/CN 

a  a  a 


(149) 


The  effect  of  a  boattail  is  to  reduce  the  normal  force  and  to  shift  the 

total  body  x  forward  (Reference  5).  The  degree  of  effect  is  a  function  of  the 
cp 

same  parameters  which  were  mentioned  in  regard  to  the  drag.  The  general  behavior 
of  the  boattail  normal  force  is  illustrated  in  Figure  63  (Reference  39).  For 


107 


SUBSONIC  -I-  SUPERSONIC 


VlMl-M 

Figure  63.  Decrease  in  Due  to  Boattail 


supersonic  flow,  the  trend  is  more  easily  seen  in  Figure  64.  For  a  given  Mach 
number  and  base  diameter,  increasing  the  boattail  length  results  in  increasing 
normal  force  losses.  Increasing  the  base  diameter  or  increasing  the  Mach  number 


has  the  opposite  effect,  resulting  in  a  restoration  of  the  normal  force.  Note 
again  that  the  Aeroprediction  Code  varies  from  the  slender  body  theory  value.  In 
this  case,  the  boattail  effect  is  underestimated  by  roughly  15  percent.  Neither 
theory  has  been  verified  by  a  detailed  comparison  with  experimental  data.  Thus, 
the  accuracy  of  these  theories  could  not  be  determi? ed.  The  location  of  the 
center  of  pressure  on  a  boattail  '.s  presented  in  Figure  65.  Again,  increasing 
the  boattail  length  for  a  giver,  oase  diameter  and  Mach  number  reduces  the  value 


of  x  /Ld  even  though  the 
cp  B 

shoulder  is  moving  rearward. 


x  location  relative  to  the  afterbody-boattail 
cp 

If  the  nose-afterbody  length  is  held  constant,  the 
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increasing  boattail  length  therefore  causes  a  strong  destabilizing  trend.  How¬ 
ever,  if  the  total  body  length  is  held  constant,  this  destabilizing  effect  is 
greatly  reduced  and  in  some  cases  eliminated. 

When  designing  a  boattail,  one  must  avoid  boattail  angles  sufficiently  large 
to  cause  flow  separation  at  the  shoulder.  This  is  particularly  the  case  for 
conical  boattails,  whereas  for  ogival  shapes  a  slightly  greater  boattail  angle 
(as  determined  by  the  local  slope)  can  sometimes  be  tolerated.  The  angle  at 
which  flow  separation  occurs  can  vary  considerably,  again  depending  on  the  local 
flow  conditions.  For  projectiles  and  other  short  bodies,  separation  appears 
generally  to  occur  for  shoulder  angles  of  6°  to  10°  or  greater.  However,  on 
boattails  following  long  afterbodies,  the  flow  has  been  known  to  remain  attached 
for  shoulder  angles  up  to  15°. 

For  spinning  bodies  such  as  projectiles,  one  should  also  examine  the  effect 
of  boattail  design  on  the  Magnus  force  and  Magnus  moment.  More  is  said  concern¬ 
ing  these  effects  in  the  following  paragraphs.  For  now,  suffice  it  to  say 
that  flight  tests  have  shown  that  boattail  design  is  important  as  far  as  the 
Magnus  force  and  moment  are  concerned,  particularly  in  transonic  flow. 

The  effects  of  nose  bluntness  on  the  normal  force  and  center  of  pressure 
are  shown  in  Figures  66  and  67  for  spherically  capped  and  truncated  conical 
noses.  The  data  presented  are  for  truncated  7.125°  semi-vertex  angle  cones  with 
a  fineness  ratio  of  4.0  calibers  for  no  blunting  (R^/Rg  =  0)  and  10°  spherically 
capped  cones.  When  adding  the  spherical  cap  or  truncating  the  nose,  the  nose 
length  was  reduced  a3  previously  described  in  the  drag  paragraph.  Overall, 
blunting  the  nose  up  to  a  bluntness  ratio  of  Rj^ref  =  • 1  has  a  negligible  effect 
on  (see  Figure  68)  and  a  small  effect  on  the  xc  location.  Theoretically, 
the  x  location  experiences  an  abrupt  shift  forward  with  initial  blunting  due  to 
truncation  (as  shown  in  Figure  69)  which  does  not  occur  when  the  spherical  cap  is 
added.  This  phenomenon  has  not  been  verified  experimentally. 

At  a  >  0,  there  is  a  nonlinear  contribution  to  the  normal  force  due  to  the 
viscous  crossflow  velocity.  The  correction  for  angle-of-attack  as  given  in 
Reference  18  [and  using  Equations  (25)  and  (26)  from  Section  2]  is 
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10°  CONE  •  TRUNCATED  TO  SHORTER  LENGTH 


3.0 


2.5 


O  R/Rf„  -  0.0 
O  R/Rr„  -  0.1 
0  R/Rf<(  -  0.2 
O  R/Rt„  -  0.4 


NOTE  -  WEIGHTED  FIT  SHOWN  BETWEEN 

SMALL  PERTUBATION  THEORY  <M„  <  2.0) 
AND  2nd  ORDER  SHOCK  EXPANSION 
THEORY  (M„  »  3.0)  (-  •  -) 


- - 


1.0  -I - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - F 
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Figure  68.  Normal  Force  for  Blunt  7.125°  (4  cal)  Cones 
Truncated  to  Shorter  Length 


Figure  69.  Center  of  Pressure  for  Blunt  7.125°  (4  cal)  Cones 
Truncated  to  Shorter  Length 


(25) 


ACM  =  C .  (M  sin  a)  S  sin2a 

N  7t  dc  c»  p 

and 

AC  =  -  ~  r)  C.  x  sin2a  (26) 

m  n  1  dc  p 


where  r|  and  are  given  in  Figure  70, 
c 

Sp  is  the  planform  area  of  the  body  in  the  plane  of  the  body's  longitudinal 
axis,  and 

S  r  is  the  reference  area, 
ref 

The  center  of  pressure  for  the  entire  body  is  then 


C  +  AC 
_  _  jn _ m 

{CP  “  "  CN  +  ACn 


(150) 


Equations  25,  26  and  145  do  not  account  f^r  any  crossflow  Reynolds  number  de¬ 
pendence.  A  more  detailed  presentation  on  the  effects  of  incidence  is  given  in 
Reference  18. 


4.1.3  DYNAMIC  AERODYNAMICS 


The  Roll  Damping  Coefficient  for  an  axisymmetric  body  is  the  result  of  the 
skin  friction  force  tangential  to  the  rolling  axis.  As  indicated  in  Section  2, 
this  coefficient  hap  been  accurately  predicted  using  the  known  values  for  a  given 
configuration  and  adjusting  these  values  through  the  relation 


(151) 


The  data  presented  in  Section  3  were  thus  obtained  and  can  be  similarly  used  if 
desired  to  estimate  the  C^  values  for  additional  configurations  so  long  as  the 
comparison  design  is  approximately  the  same  length. 
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Figure  70.  Drag  Proportionality  Factor  and 
Crossflow  Drag  Coefficient 


The  Magnus  force  and  Magnus  moment  ar>;  the  results  of  boundary  layer  distor¬ 
tion  on  a  spinning  body  at  angle-of-attack.  Usually,  empirical  methods  are  used 
to  predict  the  magnitude  of  the  Magnus  coefficients.  The  Aeroprediction  Code 
does  not  take  into  account  any  variations  in  the  nose  shape  or  length,  only  the 
total  body  length.  It  is  assumed  that  the  Magnus  moment  is  also  independent  of 
the  boattail  slope.  Sample  variations  of  the  Magnus  moment  for  bodies  with  and 
without  boattails  are  shown  in  Figures  71  and  72  which  indicate  that  increasing 
the  body  length  and/or  adding  a  boattail  increases  the  Magnus  moment.  The 
Magnus  moment  is  not  determined  for  bodies  with  lifting  surfaces.  For  missiles, 
this  is  not  significant.  For  bombs,  fin  Magnus  is  important. 

The  pitch  damping  coefficient,  C  +  C^,  indicates  the  pitching  moment 
produced  by  rotational  motion  about  the  pitching  axis  and  can  be  determined 
theoretically  as  shown  in  Section  2  (Equations  90-96).  Examples  are  given  in 
Section  3.  The  reader  is  referred  to  these  sections  for  more  information  on 
these  methods. 

The  theory  found  in  Reference  25  indicates  that  the  pitch  damping  coeffi¬ 
cient  is  strongly  influenced  by  the  nose  length  where  the  afterbody  length  is 
less  than  1.0  caliber  in  length.  For  afterbody  lengths  greater  than  1.0  caliber, 
the  main  parameter  is  the  total  body  length.  The  influence  of  the  boattail  is 
negligible  as  long  as  the  length  of  the  boattail  is  included  in  the  total  body 
length  when  determining  C  +  C^.  In  general,  increasing  the  nose  length  in 
conjunction  with  very  small  afterbodies  or  increasing  the  total  body  length  in¬ 
creases  C  .  The  adaptation  of  the  slender  body  theory  in  Reference  22  also 
indicates  that  the  increase  in  C  +  C  •  is  proportional  to  the  increase  in  the 
magnitude  of  the  pitching  moment,  C  ,  and  that,  for  configurations  with  very 

Of 

long  afterbodies  (L.  ■*  ») ,  the  value  of  C  +  C  •  approaches  the  value  of  -C  . 

°  A  ’  mq  ma  vv  ma 

Figure  73  illustrates  the  variation  in  C  +  with  body  length  as  predicted 
by  the  method  discussed  in  Reference  25. 


f 
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MACH  NUMBER 


Figure  71.  Variation  of  Magnus  Moment  for  Bodies  Without  Boattails 


Figure  72.  Magnus  Moment  Variation  for  Bodies  with  Boattails  M^  >  2.0 


4.2  LIFTING  SURFACES 


A  large  number  of  wing  alone  cases  (two  fins  together  with  total  span,  b,  no 
body  present)  were  considered  using  the  Aeroprediction  Code.  The  majority  of 
cases  represented  rather  simple  planform  shapes;  i.e.,  diamond  or  biconvex  air¬ 
foil  sections  with  a  constant  shape  and  thickness  ratio  along  the  span  and  sharp 
leading  and  trailing  edges.  Exact  linear  theory  is  applicable  in  these  cases. 
The  results  of  the  numerical  solutions  obtained  by  the  Code  are  compared  with  the 
analytical  predictions  of  the  linear  theory.  The  numerical  method  used  in  the 
Code  can  be  applied  over  a  wide  range  of  airfoil  shapes  and  plani'orms.  Modified 
double  wedge  and  biconvex  airfoils  with  blunt  leading  and  trailing  edges  and 
tapered,  swept  fins  with  variable  (but  similar)  airfoil  sections  and  variable 
thickness  with  span  are  also  possible.  However,  no  camber  or  twist  is  con¬ 
sidered.  Because  of  the  large  number  of  variations  in  geometry  possible,  the 
examples  which  follow  are  generally  limited  to  the  simple  planform  shapes.  The 
reader  is  referred  to  Section  1  of  this  report  and  to  Reference  2  for  more  de¬ 
tails  on  the  theory  utilized  in  the  Code. 

4.2.1  DRAG 

The  drag  of  a  single  fin  is  generally  small  compared  to  that  of  the  total 
configuration.  As  the  number  of  fins  increases,  the  drag  contribution  of  the 
lifting  surfaces  can  become  quite  significant,  although  rarely  dominant.  Geome¬ 
trical  factors  which  are  predominant  in  affecting  the  drag  are  aspect  ratio, 
thickness  ratio,  bluntness  of  leading  and  trailing  edges,  leading  edge  sweep 
angle,  and  airfoil  shape. 

Reference  2  contains  several  examples  which  compare  the  numerical  results 
for  the  two-dimensional  pressure  coefficient  with  that  of  the  exact  linear 
theory.  The  comparisons  verify  the  accuracy  of  the  numerical  method  for  both 
the  cnordwise  and  spanwise  distributions. 

Figures  74  and  75  show  the  wave  drag  coefficients  for  rectangular  and  delta 
planform  fins  with  sharp  LE  and  TE  for  transonic  speeds.  The  wave  drag  is  as¬ 
sumed  to  decay  linearly  from  =  1.05  to  zero  at  =  0.85.  In  general,  for 
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Transonic  Wave  Drag  Coeffici 


M»-1  M»-1 

IM«(2»/C)I»'»  |M>(2t/CII"> 


Figure  75.  Wave  Drag  Coefficient  of  a  Double-Wedge  Delta  and  Rectangular 

Biconvex  Fins  at  Transonic  Speeds 

constant  thickness  with  span,  which  is  the  case  here,  and  zero  sweep  angle,  the 
wave  drag  coefficient  based  on  the  planform  area  does  not  change  with  increasing 
span.  The  area  which  is  an  exception  to  this  rule  is  1.0  <  <  1.2.  T.n  this 
flow  region,  the  effect  of  increasing  the  aspect  ratio  is  to  increase  the  wave 
drag  coefficient.  Thickness  effects  are  also  shown  to  be  present. 

Figures  76  and  77  compare  the  supersonic  wave  drag  predictions  of  the  two 
theories  for  modified  double-wedge  and  biconvex-planform  fins,  respectively,  with 
various  sweep  angles,  y,  and  constant  chord  length.  The  agreement  is  very  favor¬ 
able.  It  can  be  seen  that  the  wave  drag  increases  with  thickness  but  decreases 
with  an  increasing  sweep  angle.  Figures  76  and  77  apply  for  fins  of  moderate 
aspect  ratios  and  moderate  supersonic  Mach  numbers.  For  higher  values  of  these 
parameters  for  rectangular  fins  with  diamond  airfoil  sections,  one  can  use  Fig¬ 
ure  78,  which  is  an  extension  of  the  curve  for  AR  tan  p  0  in  Figure  76.  One 
should  note  that  the  drag  can  change  significantly  with  sweep  depending  on 
whether  the  leading  edge  is  subsonic  (lies  forward  of  the  Mach  wave  line)  or 
supersonic  (lies  entirely  rearward  of  the  Mach  wave  line).  Reference  5  shows. 
supersonic  wave  drag  predictions  using  the  linear  theory  for  fins  with  various 
taper  ratios.  Also  shown  in  Reference  5  are  form  factors  which  can  be  applied  to 
a  variety  of  airfoil  shapes  to  adjust  the  drag  prediction  of  the  diamond  airfoil 


Figure  / 


Figure  78.  Supersonic  Wave  Drag  of  Wing  with  Double-Wedge 
Airfoil  Section,  Ct/Cr  =  1.0 

to  obtain  the  drag  for  a  fin  with  a  similar  aspect  ratio,  chord  length,  and  Mach 
number.  As  one  would  expect,  in  general,  the  greater  the  initial  slope  on  the 
airfoil  for  a  given  thickness  ratio,  the  greater  the  wave  drag  contribution. 

The  base  pressure  of  the  fin,  referred  to  as  the  trailing  edge  separation 
drag,  is  similar  to  the  three-dimensional  phenomenon  which  occurs  on  the  body. 
The  pressure  in  the. rear  of  the  fin  is  that  of  a  two-dimensional  rearward  facing 
step.  Experimental  results  (Reference  40)  for  fins  with  no  TE  slope  are  shown  in 
Figure  79.  The  data  have  been  extrapolated  for  Mw  <  1.1  based  on  the  general 
shape  of  the  three-dimensional  base  pressure  curve.  The  base  pressure  on  one  fin 
can  be  adjusted  for  TE  taper  (varying  bluntness)  by  Equation  (44)  in  Section  2. 

The  position  of  the  fin  with  respect  to  the  base  on  the  body  can  also  affect 
the  base  pressure  (of  the  body).  Also,  the  geometric  fin  parameters  which  de¬ 
termine  the  extent  of  this  effect  are  the  aspect  ratio,  thickness  ratio,  profile 
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Figure  79.  Two-dimans ional  Base  Pressure  Coefficient 


or  shape ,  sweepback  angle,  and  number  of  fins.  The  effect  of  thickness  is  shown  j 

in  Figure  80  which  indicates  that  the  effect  increases  with  an  increasing  thick-  \ 

ness  ratio  for  a  given  Mach  number.  The  other  geometric  parameters  were  not  ^ 

examined.  The  distance  or  location  of  the  fin  upstream  from  the  base  at  which  ^ 

the  fins  no  longer  affect  the  base  pressure  varies  linearly  with  t/C,  as  shown  in  3 

Figure  81  (Reference  41).  The  change  in  base  pressure  is  assumed  to  vary  linear- 
ly  from  a  maximum  when  the  TE  of  the  fin  is  flush  with  the  base  to  zero  at  the  j 

distance  indicated.  Fins  extending  downstream  of  the  base  are  assumed  to  have  ] 

the  same  effect  as  those  that  are  flush  with  the  base.  * 


The  LE  bluntness  effect,  on  the  drag  was  also  examined.  The  trends  for 
several  Mach  numbers  are  shown  in  Figure  82  for  a  biconvex  airfoil.  Bluntness  is 


effect  at  the  lower  supersonic  Mach  numbers  is  shown  to  indicate  a  stronger  de¬ 
pendence  on  the  LE  radius.  The  exact  linear  theory  does  not  hold  in  these  cases. 

The  reader  is  referred  to  Section  2  for  a  discussion  of  the  skin  friction 
drag  component. 

4.2.2  LIFT 

For  subsonic  speeds,  the  only  contributon  to  the  lift  or  normal  force  is 

assumed  to  be  due  to  the  angle-of-attack.  The  procedure  in  this  instance  (as 

shown  in  Section  2)  is  to  determine  the  normal  force  and  x  for  a  similar  wing 

cp 

of  aspect  ratio  PAR  at  Mach  number  zero  and  angle-of-attack,  pa.  The  variation 
of  for  a  diamond  airfoil  in  subsonic  flow  is  given  in  Figure  83.  For  tran¬ 
sonic  ¥low,  the  relation 

/ q  \  -  _  (1-ac)  27tAR _  (152J 

\  Na)fb  2  +  [AR2  (p2  +  tan2y^)  +  4]** 

is  used  to  compute  the  lift  curve  slope  at  the  force  break  Mach  number  and  at 

intermediate  Mach  numbers,  Ma  and  (defined  in  Section  2).  The  force  break 

Mach  number  is  assumed  to  represent  the  point  at  which  the  flow  transitions  from 

subsonic  to  transonic  flow.  The  force  break  Mach  number  varies  with  aspect  ratio 

and  thickness  as  shown  in  Figure  84-A.  The  factor,  ac,  is  a  correction  factor  to 

the  lifting  line  theory  to  account  for  the  decrease  in  the  lift  curve  slope  for 

thick  fins  as  a  function  of  thickness.  The  variation  of  this  factor  is  shown  in 

Figure  84-B  at  M^  =  Mfl  and  Figure  84-C  at  M^  =  M^.  The  further  determination  of 

in  transonic  flow  is  outlined  in  Section  2. 
a 

Figures  85  through  94  show  the  variation  in  and  x  as  a  function  of  AR, 

t/C,  and  Mach  number  for  a  rectangular  fin  with  sharp  LE  and  TE  for  transonic  and 

supersonic  flows.  Although  shown  as  a  function  of  the  transonic  similarity 

parameters,  these  curves  were  extended  to  show  the  decreasing  effect  of  the 

aspect  ratio.  For  the  rectangular  fin,  the  effect  of  AR  is  not  considered  for 

high  supersonic  Mach  numbers  where  strip  theory  is  used,  although  some  small 

effect  is  noticed  in  the  x  location.  The  effect  of  fin  thickness  in  this 

cp 

region  is  thus  more  clearly  seen.  For  M^  >  3.0,  increased  thickness  is  shown  to 
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THICKNESS  RATIO  (%  CHORD) 

Figure  84-C.  Chart  for  Determining  Lift-Curve  Slope 
Correction  Factor  at 

Figure  84.  Chart  for  Determining  Transonic  Lift-Curve 
Slope  at  Ma  and 
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Figure  85.  Normal-Force  Coefficient  Gradient  for  Rectangular  Fin 
(Diamond  Shape)  (t/C  =  .025) 


Figure  86.  Normal  Force  Coefficient  Gradient  for  Rectangular  Fin 

(Diamond  Shape)  (t/C  =  .05) 


Figure  87.  Normal-Force  Coefficient  Gradient  for  Red 

(Diamond  Shape)  (t/C  =  .1) 


Figure  90.  Center  of  Pressure  for  Rectangular  Fin 
(Diamond  Shape)  (t/C  =  .025) 


result  in  increased  C„  and  a  forward  shift  in  the  x  location.  For  low  super- 

Na  cp  y 

sonic  Mach  numbers,  the  effects  of  AK  and  yrv  are  assumed  to  dominate  the  varia- 

LL 

tion  in  C.r  and  x  and  thus  the  effect  of  thickness  on  lift  is  not  considered. 
Nfl  cp 

The  transonic  predictions  compare  well  with  the  experimental  data  correlations 
given  in  Reference  5. 

Variations  in  the  CM  and  x  with  LE  sweep  angle,  yt~,  and  taper  are  shown 

Cp  UCt 

in  Figures  95  and  96,  respectively.  Generally,  the  lift  curve  slope  effect  is 

more  noticeable  in  the  low  supersonic  flow  region  for  the  subsonic  LE  case 

(^Lg  <  M)  than  when  the  LE  is  supersonic  (Yjj.  >  p) .  (The  cases  shown  here  are 

for  unswept  TE.  The  Code  can  also  be  used  for  cases  with  TE  sweep  but,  as  these 

trends  are  similar,  the  reader  is  referred  to  P  ierence  5  for  more  detailed 

charts  produced  using  exact  linear  theory.)  Generally,  the  lift  curve  slope  is 

reduced  with  increasing  yee  a  given  taper  ratio,  \,  (increasing  AR) ,  and 

decreased  with  increasing  \  for  a  given  YEE  (decreasing  AR) .  In  supersonic  flow, 

changes  in  yee  for  a  given  \  and  Mach  number  have  very  little  effect  on  the  x^ 

location  as  a  function  of  the  root  chord,  while  increasing  \  results  in  a  forward 

movement  in  x  .  The  effect  of  thickness  was  discussed  in  the  rectangular  fin 
cp 

case  and  that  discussion  also  applies  here. 

For  complete  configurations,  the  presence  of  a  body  with  the  fin  causes 
interference  effects  on  both  the  lift  and  the  center  of  pressure.  The  inter¬ 
ference  lift  is  composed  of  three  parts:  the  lift  of  the  fin  due  to  the  presence 
of  the  body,  the  lift  of  the  body  due  to  the  presence  of  the  fin,  and  the  vortex 
lift  on  the  tail  due  to  vortices  shed  by  a  wing  or  other  forward-located  lifting 
surface.  The  methods  to  determine  these  effects  are  discussed  in  detail  in  Sec¬ 
tion  2.  The  equations  for  summarizing  the  coefficients  for  the  individual  body 
components  to  obtain  the  coefficients  for  the  total  configuration  including  the 
interference  effects  are  also  g  ven  in  Section  2. 
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Parameter  in  equation  for  akin  friction 
Aspect  ratio  (b2/S.p) 

Speed  of  sound  (ft/sec) 

Missile  acceleration  during  maneuver  (g's) 

Target  acceleration  during  maneuver  (g's) 

Parameter  in  equation  for  skin  friction,  lift  interference 
parameter 

Lift  interference  parameter,  infinite  afterbody 
Lift  interference  parameter,  no  afterbody 
Lift  interference  parameter,  short  afterbody 
Span  of  two  fins  (excluding  body)  (ft) 

Local  chord  length  (ft) 

Mean  aerodynamic  chord  (ft) 

Axial  force  coefficient 
Base  axial  force  coefficient 

Increment  in  base  force  coefficient  due  to  angle-of-attack 

Forebody  axial  force  coefficient 

Axial  force  due  to  skin  friction 

Pressure  (or  wave)  axial  force  coefficient 

Nose  separation  drag  coefficient 

Axial  force  coefficient  due  to  presence  of  a  rotating  band 

Wing  axial  force  coefficient 

Drag  coefficient  in  wind  axes 

Zero  angie-of -attack  drag  coefficient 


A-l 


Cdc 

Infinite  cylinder  crossflow  drag  coefficient 

cF,  cf 

Friction  coefficient 

CFC 

Friction  coefficient  at  critical  Reynolds  number 

C£ 

Roll  moment  coefficient 

C£p 

Roll  damping  coefficient  (9C^/3(pD/2V) ]  (rad  *) 

Cm 

Pitching  moment  coefficient 

c  +  c  • 

mq  ma 

Pitch  damping  coefficient  jo/q^SD)  {  [9My/9(qD/2V) ]  + 

[9My/9(aD/2V)]}|  (rad-1) 

C 

ma 

Pitching  moment  coefficient  derivative 

C 

m6 

Pitching  moment  coefficient  derivative  due  to  fin  deflec¬ 
tion, 

CN 

Normal  force  coefficient 

CNB 

Body-alone  normal  force  coefficient 

CN 

NB(F) 

Normal  force  coefficient  correction  for  interference  due  to 
body  in  presence  of  a  fin 

NF(B) 

Normal  force  coefficient,  correction  for  interference  due  to 
fin  in  presence  of  a  body 

CNTi 

th 

Isolated  1  fin  (in  presence  of  the  body)  normal  force 

coefficient  (cruciform  configuration) 

CNT(V) 

Tail  normal  force  coefficient  correction  for  canard  shed 
vortex 

^N (TOT) '  CN(B+T) 

Body-tail  total  normal  force  coefficient  at  small  a  and  high 
a,  respectively 

CN 

a 

Normal  force  coefficient  derivative 

C 

V 

Magnus  moment  coefficient  derivative  { (1/q^SD)  [92Mz/9(pD/2V)dof] 
(rad  *) 

C 

P 

Pressure  coefficient 

A-2 


Ol 


CpL’  CpU 


Pressure  coefficients  for  lower  and  upper  surfaces  of  the 
wing,  respectively 

Base  pressure  coefficient 


Fin  root  chord  (ft) 

Distance  from  wing  root  leading  edge  to  first  surface  dis¬ 
continuity  downstr  'in,  parallel  to  freestream  (ft) 

Distance  from  wing  'ot  trailing  edge  to  first  surface  dis¬ 
continuity  upstream  parallel  to  freestream  (ft) 

Fin  tip  chord  (ft) 

Distance  from  wing  tip  leading  edge  to  first  surface  dis¬ 
continuity  downstream,  parallel  to  freestream  (ft) 

Distance  from  wing  tip  trailing  edge  to  first  surface  dis¬ 
continuity  upstream,  parallel  to  freestream  (ft) 

Parameter  in  equations  for  skin  friction,  transonic  nose 
normal  force 


C2 

CY 

CG 

D 


Parameter  in  equations  for  skin  friction,  transonic  nose 
normal  fo-ce 

Constant  in  determining  pitch  damping  coefficient 

Center  of  gravity 

Body  reference  diameter  (ft) 

Mean  body  diameter 

Mean  body  diameter  near  a  canard  and  tail,  respectively 

Parameter  in  equation  for  skin  friction 

Parameter  in  equation  for  skin  friction 

Mach  wave  angle  proportionality  factor 

Spanwise  location  of  panel  vortex 

Dynamic  pressure  ratio  across  a  bow  shock 


A-3 


i 


Acceleration  due  to  gravity  (ft/sec2) 

Rotating  band  height  above  body  (cal) 

Initial  heading  error  (rad) 

Body-tail  normal  force  coefficient  interference 
Specific  impulse  (sec) 

Body  polar  inertia  (slug-ft2) 

Body  transverse  (pitch  axis)  inertia  (slug-ft2) 
Tail  interference  factor 
Hypersonic  similarity  parameter 


KB’  KC’  KB(F)’ 
KF(B)’  kB(F)’  ^(B) 


LA’  LB’  LN 


LB(F)’  LF(B) 


M  .  M, 
a  ’  b 


Interference  parameters 


Body  length  (cal) 

Lengths  of  afterbody,  boattail,  and  nose,  respectively  (cal) 
Body  alone  lift 
Interference  lift 

Combination  lift  of  wing  plus  body 
Fin  alone  lift 
Comparable  body  length 
Mach  number 

Intermediate  transonic  Mach  numbers 
Force  break  Mach  number 

Mach  number  separating  low  supersonic  and  high  supersonic 
inviscid  body  alone  and  wing  wave  drag  computations 

» 

Mach  number  separating  low  supersonic  and  high  supersonic 
inviscid  lift  and  dynamic  derivative  computations 


M* 


m 


mHE 

N 

P 

1 

P 

R 

"be 

"n 

"nc 

rne 

"nl 


"x 


Mach  number  related  to  nose  length  used  for  interpolations 
Mass 

Miss  distance  due  to  heading  error  fp'eters) 

Effective  navigation  ratio 
Roll  rate  (rad/sec) 

Pitch  rate  (rad/sec) 

Dynamic  pressure 
Body  radius  (cal) 

Mean  local  body  radius  (cal) 

Radome  boresight  error 

Body,  nose  spherical  radius  (cal) 

Critical  Reynolds  number 

Effective  radius  of  spherical  cap  (cal) 

Reynolds  number  based  on  body  length 
Body  radius  at  base  of  spherical  cap  (cal) 

Body  s  ‘ace  slope  at  x=0 

Body  surface  slope  at  nose-afterbody  junction 
Range  (ft) 

Area  (unit2) 

Body  surface  area  (cal2) 

Dynamic  stability  factor 

Gyroscopic  stability  factor 

Body  longitudinal  planform  area  (cal2) 

Total  planform  area  (ft2) 


W 


w 


s 

st 

T 

T 

R 

T 

w 

t 

t 

t„ 


V 

V 


m 


W 

WT 


Planform  area  for  two  tail  fins  (ft2) 

Wetted  area  (ft2) 

Distance  along  linear  segment 

Point  on  body  at  which  perturbation  solution  begins 
Origin  of  source  or  doublet  distribution 
Temperature  (°R) 

Thrust  (lb) 

Wall  temperature  (°R) 

Local  thickness  of  lifting  surface  (ft) 

Body  surface  match  point  for  spherical  cap  (cal) 

Time  of  flight  (sec) 

Total  missile  maneuver  time  constant 
Lifting  surface  thickness  at  root  (ft) 

Lifting  surface  thickness  at  tip  (ft) 

Velocity  (ft/sec) 

Closing  velocity  (ft/sec) 

Missile  velocity  (ft/sec) 

Initial  velocity  (ft/sec) 

Weight  (lb) 

Launch  weight  (Lb) 

Perturbation  velocity  (ft/sec) 

Stretched  coordinate  at  which  singularity  exists 

Body  coordinate  (parallel  to  body  axis)  or  lifting  surface 
chordwise  coordinate  for  fin  only  (cal) 

Body  coordinate  including  spherical  cap  (cal) 


A-6 


AFT 


eg 


X  ,  X 

c’  w 


cp 

xcpB 

^cp^Ti 

^cp^TOT 


cpi 


te 


y 

y  . 

Jcpi 

a 

P,  P 

r 

Y 

Y 

*G 

Yl >  Y2»  Y3>  Y4 


Distance  from  tail  trailing  edge  to  base  of  body  (cal) 

Distance  from  nose  to  moment  reference  (cal) 

Distances  to  fin  apex  for  canard  and  tail,  respectively, 
from  nose  (cal) 

Distance  from  moment  reference  to  center  of  pressure  (cal) 
Body-alone  x  from  nose  (cal) 

Isolated  i*"*1  fin  (cruciform  configuration)  x  from  nose  (cal) 

Total  body-tail  x^  from  nose  (cal) 

Interference  x  from  nose  (cal) 
cp 

Body  planform  area  center  location  (cal) 

Distance  defined  in  Equation  (125) 

Stretched  coordinate  at  which  singularity  exists,  functional 
parameter,  radius  coordinate 

Spanwise  distance  on  lifting  surface  (ft) 

Spanwise  location  of  center  of  pressure  for  isolated  i^  fin 
(ft) 

Angle-of-attack  (rad) 

Compressibility  similarity  parameter,  -  1  | 

Second  derivative  of  C  with  respect  to  a  evaluated  at 
a=0 ,  6=0  p 

Calorically  perfect  gas  heat  capacity  ratio 

Flight  path  angle  of  missile  rate-of-change  (deg/sec) 

Streamwise  circulation  per  unit  length 

Sweep  angles  for  local  surface  discontinuities  on  a  fin, 
functional  parameter  (deg,  rad) 

Leading  edge  sweep  angle  of  lifting  surface  (=Yi)  (deg,  rad) 

Incremental  value,  second  derivative  of  C  with  respect  to 
a  evaluated  at  a=0,  0=n/2  P 


A-7 


6 

fi* 

6F 

£ 


n 

e 


e 


B 


0 


c 


0* 


A 


A 


M 

i 

P 


Jit 


it 


tr 


tt 


Local  body  angle  (rad) 

Nose  slope  angle  (rad) 

Fin  control  deflection  angle  (deg) 

Angle  associated  with  body  slope  (rad) 

Angle  of  nose-afterbody  junction  (rad) 

Initial  body  surface  angle  following  spherical  cap  (rad) 
Crossflow  drag  proportionality  factor,  functional  parameter 
Azimuthal  angle  from  the  leeside,  roll  angle  (rad) 

Boattail  or  flare  angle  (deg) 

Cone  half  angle  (deg) 

Equivalent  cone  angle  for  nose  of  given  length  (rad) 

Parameter  proportional  to  derivative  of  C  with  respect  to 
a  for  a=0,  0=0,  functional  variable  p 

C  /C 
V  r 

Mach  wave  angle,  sin 
Functional  parameter 
Density  (slug/ft3) 

Root  leading  edge  radius  (cal) 

Tip  leading  edge  radius  (cal) 

Root  trailing  edge  radius  (cal) 

Tip  trailing  edge  radius  (cal) 

Area  of  influence  on  lifting  surface,  stress  (psi) 

Time  (sec) 

Airframe  time  constant 
Guidance  filter  time  lag  constant 


A-8 


Total  missile  guidance  time  constant 
Tracking  loop  time  constant 
Disturbance  potential 
Functional  parameter 


Differentiation  with  respect  to  x 
Span  including  body 

Afterbody 
Asymptotic  value 
Base 
Body 

Conical  value 

Local 

Nose 

Ogival  value 
Pressure 
Rotating  band 
Reference  value 
Standard  value 
Tail  value 
Wing 


Wave  drag  component 


Cont'd) 


Stagnation  value  (when  referring  to  pressure,  except  on 
pages  24-25  where  it  refers  to  the  M  =0  case),  value  at  a=0 
(drag  only)  06 

Values  upstream  and  downstream  from  a  body  corner, 
respectively 

Partial  differential  with  respect  to  a  (a+0) 

Freestream  condition 

Parameter  in  determination  of  C  +  C  •  using  LMSC  method 
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